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Abstract 

Consider a system of N bosons in three dimensions interacting via a repulsive short range 
pair potential N 2 V(N(xi — Xj)), where x = (x%, . . . ,xn) denotes the positions of the particles. 
Let Hn denote the Hamiltonian of the system and let tpN,t be the solution to the Schrodinger 
equation. Suppose that the initial data tpN,o satisfies the energy condition 

{ipN,o, H N ifj Nfi ) < CN . 

and that the one-particle density matrix converges to a projection as N — > oo. Then, we prove 
that the fc-particle density matrices of ipN,t factorize in the limit N — > oo. Moreover, the one 
particle orbital wave function solves the time-dependent Gross-Pitaevskii equation, a cubic non- 
linear Schrodinger equation with the coupling constant proportional to the scattering length of 
the potential V. In [12], we proved the same statement under the condition that the interaction 
potential V is sufficiently small; in the present work we develop a new approach that requires no 
restriction on the size of the potential. 

1 Introduction 

We consider a bosonic system of N particles with a repulsive interaction. The states of the system are 
given by elements of the Hilbert space L 2 (IR 3Ar ), the subspace of L 2 (M. 3N ) consisting of permutation 
symmetric wave functions. We are interested in describing the time evolution of special initial wave 
functions tpN G L 2 (M 3Ar ) that exhibit complete Bose-Einstein condensation. 

For a given wave function tpN, we define the density matrix 77V = \iPn){iPn\ associated with ^jv 
as the orthogonal projection onto ipN- Moreover, for k = 1, . . . , N, we define the /c-particle marginal 
density 7^ , associated with ipw, by taking the partial trace of 7^ over the last (N — k) variables. 
In other words, 7^ is defined as a positive trace-class operator on L 2 (M 3fc ) with kernel given by 

7^' ) (x fc ;x' fc ) = j dx N „ k il' N (x k ,x N ^ k )ip N (x k ,x N - k ) . 
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Here and in the following we use the notation x& = (x%, . . . ,Xk), x' fc = (x[, . . . ,x' k ), and xjy^k = 
(xk+i, • • • , xn); we will also denote x = (x\, X2, ■ ■ ■ ,%n)- A sequence {V'jvlweN, with ifi^ S L;?(R 3Ar ) 
for all N, is said to exhibit complete Bose-Einstein condensation in tp € L 2 (R 3 ) if 

^ I^X^I asiV^oo (1.1) 

in the trace-norm topology (here and in the following \ip){(p\ indicates the orthogonal projection onto 
if). Physically, complete Bose-Einstein condensation means that all particles in the system, apart 
from a fraction vanishing as N — > 00, are described by the same one-particle wave function 92, known 
as the condensate wave function. Note that (jl.ip implies that 

(1.2) 

for all k > 1, as was first proven by Lieb and Seiringer in |18j . 

The time- evolution of N boson systems is governed by the Schrodinger equation 

id t ipN,t = H N ip N j (1.3) 

with the Hamiltonian operator 

N 

j=\ i<j 

Here and in the following we are going to use the convention Aj = A Xj and Vj = V x .. We consider 
the scaling introduced by Lieb, Seiringer and Yngvason in [20] for the interaction potential Vjy, i.e. 
we fix a non-negative potential V, and then we rescale it by defining V^(x) = N 2 V(Nx). This 
scaling is chosen so that the scattering length of Vm is of the order 1/N. 

We recall that the scattering length associated with a potential V decaying sufficiently fast at 
infinity (V has to be integrable at infinity) is defined through the solution of the zero-energy scattering 
equation 

(_A + M/ = (1.5) 

with boundary condition f{x) — >• 1 as \x\ — ► 00. We usually write / = 1 — w. The scattering length 
of V, which is a measure of the effective range of the interaction, is defined by 

ao = lim |a;|iy(x). 

|x|^oo 

An equivalent definition of the scattering length is given by the formula 

J dx V(x)(l - w(x)) = 8vra . (1.6) 

By these definitions, it is clear that, if ao denotes the scattering length of the potential V, then the 
scattering length of the scaled potential Vjsr is given by a = clq/N. 

Our main result is as follows; suppose that the family of wave functions {VwjiVeN exhibits com- 
plete Bose-Einstein condensation (jl.lh with some if G ^(R 3 ) and assume its energy per particle 
to be bounded (in the sense that {i^n, -ffjv^Tv) < CN for all N). Denote by ipN,t the solution of 
the Schrodinger equation (|1.3p with initial data ipN,o = ^N- Under appropriate conditions on the 
potential V, we show that, for every time t 6 R, the family {Vw^jiVeN still exhibits complete con- 
densation, and that the condensate wave function evolves according to a cubic nonlinear Schrodinger 
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equation known as the Gross-Pitaevskii equation. In other words, if denotes the one-particle 
marginal density associated with ipN,t, we prove that 

tSH^X^I (1-7) 

as N — > 00, where <pt is determined by the nonlinear Gross-Pitaevskii equation 

id t <pt = —Atpt + 8irao\ipt\ 2 (ft (1-8) 

with initial data tpo = (p. The cubic non-linear term expresses the local on-site self-interaction 
of the condensate wave function. Due to the strongly localized interaction, the many-body wave 
function develops a singular correlation structure on the scale 1/N. As N goes to infinity, this short 
scale structure produces the scattering length as a coupling constant in the limiting Gross-Pitaevskii 
equation. 

This result gives a mathematical description of the dynamics of initial data exhibiting complete 
Bose-Einstein condensation. The simplest example of such initial data are product states ipN = 
(p® N , for arbitrary ip £ L 2 (R 3 ). Physically more interesting examples of complete Bose-Einstein 
condensates are the ground states of trapped Bose gases. A trapped Bose gas in the Gross-Pitaevskii 
scaling is an TV-boson system described by the Hamiltonian 

N N 

H T P = Yl (~ a *j + ^xt(^o) + Yl v ^ - x i) . ( L9 ) 

j=l i<j 



where V ex t(x) > with V ex t(x) — > 00 as |x| — ► 00 is a confining potential. The Hamiltonian (11. 9j) 
describes therefore a system of iV particles confined by the external potential V ex t into a volume of 
order one, interacting through a potential with effective range of the order 1/N. Since the typical 
distance between neighboring particles is of order iV^ 1 / 3 , and thus much bigger than the effective 
range of the interaction, (|1.9p describes a very dilute system. 



In [20j, Lieb, Seiringer and Yngvason studied the ground state energy En of the Hamiltonian 
(jl.9p with a spherically symmetric interaction, V(x) = V(\x\), and they proved that 

lim — — = min £qp(lp), 

where <?gp denotes the so-called Gross-Pitaevskii energy functional 

Sgp(<p) = J dx(|V^| 2 + K xt (x)|^(x)| 2 + 4vrao|^(x)| 4 ) . 

In [18] , Lieb and Seiringer proved then that the ground state of (jl.9p exhibits complete Bose-Einstein 
condensation in the minimizer <Pgp of the Gross-Pitaevskii energy functional £gp- hi other words, 
they proved that, denotes the ground state vector of (jl.9p . then the corresponding one-particle 

marginal density 7^J trap satisfies 

as — ► cxd, in the trace- norm topology. A survey of results concerning the ground state properties 
of the dilute bosonic gases can be found in |19j . 

Since Vv ap describes a complete Bose-Einstein condensate, we can apply (|1.7p to study its time 
evolution with respect to the Hamiltonian (jl.4p ; it follows that, for every fixed t £ R, e- iHNt *p% ap 
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exhibits complete Bose-Einstein condensation in the one-particle state described by the solution 
(ft of the Gross-Pitaevskii equation (jl.8p with initial data </>gp- This result gives the mathematical 
description of recent experiments (initiated in [TJE]) where the dynamics of Bose-Einstein condensates 
has been observed. 

We already proved (|1.7p in [12] (some partial results were previously obtained in [E]), under the 
assumption of a sufficiently weak potential V. More precisely, in [12j we required the dimensionless 
parameter 

p= sup \x\ 2 V(x) + [ PjV(x) (1.10) 
xem. 3 J \ x \ 

to be sufficiently small. In the present paper, we remove this condition and consider arbitrary 
repulsive potentials V > 0, with the fast decay property V(x) < C(x) _CT , for some a > 5 (here 
(x) = (1 + x 2 ) 1 / 2 ). 

The removal of the smallness condition requires completely new ideas. The main challenge in 
the derivation of the Gross-Pitaevskii equation is to identify and resolve the short scale correlation 
structure in the iV-body wave function. In [TJ] we achieved this by locally factoring out the solution 
of the zero energy scattering equation fjl .5f) . This approach, however, is not sufficiently precise 
for a large interaction potential. In the present paper we propose an intrinsic characterization of 
the correlation structure in terms of the two-particle scattering wave operator. More precisely, we 
prove that the action of the wave operator in the relative coordinate gularizes ipN t m 

this variable. This regularity is essential to control the convergence of the many-body interaction 
to the local on-site nonlinearity in the limiting equation (|l,8p . An a-priori estimate leading to this 
regularity will be obtained from the conservation of the second moment of the energy, i.e. from 
(ipN,t) H%TpN t t) = (tpN, HfjipN}- This a-priori bound, however, only controls a specific combination 
of two derivatives, V Xi ■ V Xj , acting on the regularized wave function. We thus need to establish a 
new Poincare-type inequality involving this combination of derivatives only. In the next section we 
discuss the main ideas of our new approach. 



2 Resolution of the correlation structure for large potential 



As in [12j . the general strategy of our proof is based on the study of solutions of the BBGKY 

(k) 

hierarchy of equations for the marginal densities ^ Nt associated with the solution of the iV-particle 
Schrodinger equation (II. 3j) : 



(fey 
r N t 



i=i 



i<3 



(2.1) 



+ (N - k) E Tr fc+i Vn{xj - x k+ i),T N 

3=1 



(fc+i) 
t 



for k = 1, . . . , N. Here Tr^+i denotes the partial trace over the {k + l)-th variable, and we use 

(k) 

the convention that ") N t = for k = N + 1. The main observation is the fact that limit points 

{7oo\}fc>i °f the families {^\}^ = i (with respect to an appropriate topology) are solution of the 
infinite hierarchy of equations 



idtit\ = E 



(k) 

ii 7oo,t 



+ 8vra E Tr fe+i S ( x 3 ~ x k+i),loo,t 
3=1 



(fc+1) 



(2.2) 
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It is easy to check that the product ansatz 7^ = \<Pt) (<Pt\® k satisfies (|2.2p if and only if ipt solves the 
Gross-Pitaevskii equation (jl.Sp . Therefore, to conclude the proof of (|1.7p . it suffices to show that: 1) 
every limit point of the family {7jv+}fcLi is a solution of the infinite hierarchy (|2.2p and 2) the solution 
to ()2.2p is unique. This strategy has already been used to derive the non-linear Hartree equations for 
the effective dynamics of so called mean-field systems, see J27J H2J S E] , to derive the cubic non-linear 
Schrodinger equation with different (and simpler) scalings of the interaction potential, see [8| 111], 
and to derive the non- linear Schrodinger equation in a one-dimensional setting see [H [2] . We remark 
that the first derivation of the Hartree equation was obtained with a different method in |17[ HH] . 
With this method the speed of convergence was estimated recently in |25j . 

In all works based on the BBGKY hierarchy, the key step consists in finding an appropriate norm 
and space of density matrices to work with. On the one hand, the topology has to be sufficiently 
strong to guarantee the convergence of (12. ip to (|2.2p and the space has to be sufficiently small to 
guarantee the uniqueness of the solution to (|2.2p . On the other hand the norm defining this space 
cannot be too strong since we have to prove, via an a-priori bound, that limit points of the sequence 
7^j belong to this space. 

In [12] , we use an appropriate Poincare-type inequality to prove the convergence of (|2.ip to (|2.2p . 
To do that, we would need a control on a mixed Sobolev norm on j Nt and 7^ t with at least two 
derivatives (note that the commutator with the delta function in (|2.2p is even ill-defined unless some 
regularity is known on "f^\ )■ However, due to the singularity of the interaction potential, it turns 
out that the solution of the Schrodinger equation ^jv,t develops a short scale correlation structure, 
living on a length-scale 0(1/N), which causes the Sobolev norms with two or more derivatives to 
blow up as N — > 00. 

Instead of considering derivatives of ipN,t, we prove therefore an a-priori bound of the form 
/ |ViV2<^i2,AT (*)| 2 < Con the iV-body function 0i2,jv(i) = "0jv,t/(l — ^12) obtained from the original 
wave function after factoring out the singular short scale structure. Here 1 — w\2 = 1 — wn(x\ — X2), 
where /jv = 1 — wn is the zero-energy scattering solution to (—A + (1/2)V/v)/jv = 0. Note that, 
by simple scaling, wn(x) = w(Nx), where f(x) = 1 — w(x) is the zero-energy scattering solution 
to the unsealed equation (—A + (1/2)V)/ = 0. It turns out that J \V\V2^N,t\ 2 grows with N, 
but J |ViV2<^i2,at (*)| 2 remains bounded. Although ipN,t an d 4>i2,N(t) behave very differently in the 
mixed Sobolev norm, their difference in L 2 -norm vanishes in the N — ► 00 limit due to the scaling 
wn(x) = w(Nx). This allows us to obtain control on the mixed Sobolev norm of 700,4 despite the 
fact that it is defined as a limit of 7jvt only in the weak topology of trace class operators. Moreover, 
the boundedness of J |ViV2<^i2,JV (*)| 2 explains the emergence of the scattering length in (12. 2p . 

The proof of f |ViV20i2,A r (*)| 2 < C relies on the conservation of along the time evolution 
and on the key inequality 

2 



{^ Nl (H N + Nf^ N )>CN 2 



'1 - WV2 



(2.3) 



valid for all G L" 

To show the uniqueness of the solution of the infinite hierarchy (|2.2p . on the other hand, more 



information on the limiting densities 7^ were needed; more precisely, uniqueness was proven in [TT] 
under the assumption that 

Tr (1 - Ax) ... (1 - A,hS < C k (2.4) 

for all A; > 1. Because of the singular short scale structure characterizing the solution of the 
Schrodinger equation for finite N, the densities ^ N t do not satisfy (|2.4p . To circumvent this problem, 
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we derived in [12] a higher order energy estimate of the form 

(4> N ,(H N + N) k %b N )>C k N k J dxG^x) |Vi...V fc Viv(x)| 2 (2.5) 

for all k > 1, where Ofc_i(x) is a cutoff function vanishing (up to exponentially small contributions) 
in regions where a second particle comes close to xj, for some j < k — 1 (see Section [7] for the precise 
definition of the cutoff). This higher order energy estimate provides a control on the L 2 norm of 
the mixed derivatives Vi . . . V^iVt restricted on regions with no other particle close to x%, . . . , x^-i- 
Choosing @k-i to vanish in a sufficiently small volume, it was possible to remove the cutoffs in the 
weak limit N — > oo and to obtain the a-priori bounds (|2,4p on the limit points {ljxt}k>l- 

The estimates (|2.3p and (|2.5p were therefore the two main ingredients used in [12] to control the 
singularity of the interaction potential Vn and the singular short scale structure of ipN,t- In [32], both 
these estimates heavily relied on the smallness of the parameter p introduced in (jl.lOp . Therefore, 
although the general strategy of the current paper is similar to the one used in [12], the removal of 
the smallness condition on p requires completely new ideas, which we now explain. 

In [12], the energy estimate (|2.3p was obtained from an identity of the form 

{iP N , H 2 N iP N ) = N(N -I) j |ViV 2 0i 2 | 2 (l - w 12 ) 2 + N(N -I) j Vi0i2 • gi 2 V 2 <pi 2 + , (2.6) 

where c/>i 2 = — w\ 2 ). Here g\2 is an explicit matrix involving the Hessian of w\2 and squares 

of its first derivatives and Q contains irrelevant terms. The following bound is essentially optimal 
for the size of g\2 ■ 

\gi2\ < | [2 • 

\X\ — X2\ z 

The best strategy is then to estimate 

j Va<^i 2 • g 12 V 2 (f>i 2 > -Cp J r^j^p (|Vi</»i 2 | 2 + |V 2 12 | 2 ) > -Cp j |ViV 2 0i 2 | 2 , (2.7) 

where we used Hardy's inequality in the last step. This term can be absorbed in the first (positive) 
term in (12.61) only if p is sufficiently small. A closer inspection of the structure of the terms in Q 
reveals that, although some of them are positive, they cannot compensate for the negative term 
(|2.7p . After estimating them from below at the expense of adding lower order terms to Hfj, we get 
the desired bound 

J |ViV 2< /» 12 | 2 < CN- 2 (ip N , (H N + N) 2 ^ N ) . (2.8) 

The first main idea of the current paper is to prove the following replacement for (|2.8p (see 
Proposition 15.21) 

' | Vi • V^Vwl 2 < CN~ 2 (ip N , (H N + N) 2 ^ N ) , (2.9) 



where W\2 denotes the wave operator of the one particle Hamiltonian — A + \Vn acting on the 
difference variable x\ — x 2 . Note two main differences between (|2.8[) and (|2.9j) . First, W* 2 ipN 
replaces the function 4>\ 2 ,n = "0iv/(l — ^12) that can be considered as a first order approximation 
to W* 2 ipN- Second, instead of controlling the full mixed second derivative, ViV 2 , as in (|2.8|) . the 
new bound controls only Vi • V 2 . Although the control on Vi • V 2 is in general much weaker than 
the control on ViV 2 , in the radial direction of the relative coordinate x± — x 2 the new bound is as 
strong as the former one. 
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Both differences cause substantial difficulties in proving that (|2.ip indeed converges to (|2.2p . 
First, instead of working with a relatively explicit function w\ 2 and using its fairly straighforward 
properties summarized in Lemma 5.1 of [T2], now analogous properties have to established for the 
wave operator. Second, the lack of the control on the full mixed second derivatives impedes using 
the Sobolev-type operator inequality 

V(x 1 -x 2 )<C\\V\\ 1 (l-A 1 ){l-A 2 ) 

that was crucial in controlling many error terms. We have found a replacement for this inequality 
(Lemma 110. ip : 

V( Xl - x 2 ) < C\\V\\ X ((Vi • V 2 ) 2 - A x - A 2 + l) (2.10) 

that uses only the Vi • V 2 combination in the highest order term. Similarly to Lemma 7.2 of [12], we 
are also able to improve (|2.10p to a Poincare-type inequality (ip, (h a (xi — x 2 ) — 5(xi — x 2 ))ip) — > 
as a — > 0, where h a (x) = a~ 3 h(x/a) (with J h(x)dx = 1) is an approximate Dirac delta function on 
scale a (see Lemma Il0.2p . This is necessary to control the convergence of (|2.ip to (|2.2p . 

The second main novelty of this paper is a proof for the higher order derivative estimates (12. 5p 
in the large potential regime (Proposition 17. Ip . Although the main conclusion is the same as Propo- 
sition 5.3 in |12j . the proof does not require the smallness of p. The proof in |12j started from the 
trivial energy estimate ^2f = i — A? ^ Hn and the estimate (|2.8p on Hfj and it used an induction on 
the exponent to pass from (Hn + N) k to (Hjy + N) k+2 . The step two induction allowed us to start 
the induction estimate with 

(H N + N) k+2 > (H N + N)Vl . . . V^eV fc . . . Vi(Hjv + N) , (2.11) 

and then we commuted the two Hn + N factors through to the middle and used the (-ffjv + -^O 2 
estimate (12 .8p . The weaker estimate (I2.9p however does not allow us to regain a control on full 
derivatives, so pursuing this path would, at best, yield a control on some complicated combinations 
of partial derivatives. 

In this paper we establish higher order derivative estimates by a step one induction, i.e. passing 
from (H N +N) k to (H N + N) k+1 . This eliminates using the Hfj estimate with incomplete derivatives, 
but the price is that instead of (|2.1ip we have to work with 

(H N + N) k+1 > (H N + N)V 2 V! . • • V^eV fc . . . Vi(H;v + N) l l 2 (2.12) 

in the induction step, i.e. we have to commute the square root of the Hamiltonian through the 
derivatives. The technical complications involved with the square root turn out to be reasonably 
shortly managable, so that this new method actually provides a simpler proof than in [12] for the 
higher energy estimates even when p is small. 

Notation. Throughout the paper we will use the notation x = (a?i, . . . ,xn) £ M. 3N , and, for 
k = l,...,N, x fc = ( Xl ,...,x k ), x' fc = (x[,...,x' k ) e M 3k , xjy.fc = (x k+1 ,...,x N ) € R^ N ~ k \ The 
notation || • || indicates the L 2 -norm if the argument is a function and it denotes the operator norm 
(from L 2 to L 2 ) if the argument is an operator. For 1 < p < oo, ||/|| p indicates the L p -norm of /. 
Moreover, we will use Vj and Aj as shorthand notations for and, respectively, A Xj . If A is an 
operator acting on L 2 (M 3fc ), we will denote its kernel by ^4(xfc;x' fc ). The letter C denotes universal 
constants that may depend on V and on the H 1 -norm of the initial one particle wave function ip, 
but is independent of N. 
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3 Main Theorem 



The following theorem is the main result of this paper. 

Theorem 3.1. Suppose that V > 0, with V(—x) = V{x) and V{x) < C{x)~ a , for some a > 5, and 
for all x G M 3 . Assume that the family tp^ G L^(R 3N ), with \\iPn\\ = 1 for all N, has finite energy 
per particle, in the sense that 

{^ N ,H N i> N ) <CN (3.1) 

and that it exhibits complete Bose-Einstein condensation in the sense that the one-particle marginal 
associated with ipN satisfies 

7jV — * l < / 9 )(9 9 l i n the trace norm topology as N —* oo (3-2) 
for some tp G ff 1 (M 3 ). Then, for every k > 1 and t G R, we /iave 

as N ^ oo, in the trace norm. Here tp t is the solution of the nonlinear Gross- Pitaevskii equation 

idtft = -A(£i + 8irao\pt\ 2 *Pt (3.3) 

with initial data <pt=o = P- 

Remark. Note that the condition V{x) < C(x)~ a for some a > 5 and for all x G R 3 is only 
required to apply the result of Yajima in [29], which guarantees that the wave operator W associated 
with the Hamiltonian f) = —A + (1/2)V maps L P (M 3 ) into itself, for all 1 < p < oo (see Proposi- 
tion [5J])- If one knows, by different means, that ||W||£p_>£p < oo (it suffices to know it for p = 1 
and p = oo), then it would be enough to assume that V G L 1 (M 3 , (1 + |x| 2 )dx) n L 2 (R 3 ,dx). 

Remark. Compared with our previous result in [12], we do not require here the potential V to 
be spherical symmetric (we only need that V(—x) = V{x)). 

Remark. The fact that tp G ^(R 3 ) does not need to be assumed separately, since it already 
follows from the assumption (|3.ip . 



To prove this theorem we will use an approximation argument and the following theorem, which 
proves Theorem 13.11 for a smaller class of initial iV-particle wave functions. 

Theorem 3.2. Assume the same conditions on the potential V as in Theorem \3.1[ Suppose moreover 
that \\7 a V(x)\ < C for all multi-indices a with \a\ < 2. Assume that the family ip^ G L 2 (R 3N ), with 
\\iPn\\ = 1; is such that 

(iP N ,H N ^ N ) <C k N k (3.4) 

for all k G N, and that 

In ~ * I^X^I * n the trace norm topology as N — > oo (3.5) 
for some tp G -?/ 1 (M 3 ). Then, for every k > 1 and t G R 

as N — > oo, in the trace norm. Here tp t is the solution of the nonlinear Gross-Pitaevskii equation 
\3. 3\) with initial data tpt=o = P- 
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4 Proof of the main theorem 

In this section we present the proof of Theorem 13,21 and we show how it implies Theorem 13-H making 
use of several key proposition, whose proof is deferred to subsequent sections. 

We start by defining an appropriate space of time-dependent density matrices. To use Arzela- 
Ascoli compactness argument, we will need to establish the concept of uniform continuity in this 
space, thus we have to metrize the weak* topology. 

Let /Cfc = /C(L 2 (IR 3fc )) denote the space of compact operators on L 2 (R 3fc ) equipped with the 
operator norm topology and let C\ = C 1 (L 2 (R 3k )) denote the space of trace class operators on 
L 2 (M 3fc ) equipped with the trace norm. It is well known that C\ is the dual of JCk- Since Kk is 
separable, we can fix a dense countable subset of the unit ball of we denote it by {Ji^}i>i G 

with || J- || < 1 for all i > 1. Using the operators J- we define the following metric on the space 
C\ = £ 1 (L 2 (R 3fe )): for 7(^7^) e £\ we set 

oo 

%(7 {fc) ,7 {fc) ):=E 2 " 2 

i=l 

Then the topology induced by the metric % and the weak* topology are equivalent on the unit ball 
of C\ (see [26], Theorem 3.16) and hence on any ball of finite radius as well. In other words, a 
uniformly bounded sequence 7^ G C\ converges to G C\ with respect to the weak* topology, if 
and only if ?]k(l$ il^) —> as N — > 00. 

For a fixed T > 0, let C([0,T],££) be the space of functions of t G [0, T] with values in C\ which 
are continuous with respect to the metric rjk- On C([0, T],C\) we define the metric 

Vk(l (k) (-),l {k) (•))== sup Vk(l (k) (t),^(t)). (4.2) 
te[o,T] 

Finally, we denote by r pro( i the topology on the space © j!c > 1 C([0, T], C\) given by the product of the 
topologies generated by the metrics % on C([0, T],C\). 

Proof of Theorem \3.2l The proof is divided in four steps. 

Step 1. Compactness ofT]y,t = {lN\}k>i- We fix T > and work on the interval t G [0, T]. 
Negative times can be handled analogously. 

In Theorem 16.11 we show that the sequence = {7^j}/t>i G C([0, T], C\) is compact 

with respect to the product topology r pro d defined above (we use the convention that j K N J t = if 
k > N). It also follows from Theorem l6 . 1 1 that any limit point = {7^\}jfc>i G ©fc>! C([0, T], C\) 
is such that, for every k > 1, 7^^ > 0, and 7^^ is symmetric w.r.t. permutations. 

Using higher order energy estimates from Proposition 17. 1\ we show in Theorem 17.31 that an 
arbitrary limit point = {7^}fc>i of the sequence \ (with respect to the topology r proc j) is 
such that 

Tr(l-A 1 )...(l-A fc ) 7 S<C fe (4.3) 
for every t G [0, T] and every ft > 1. 



Tr J; 



^( fc ) 



(4.1) 
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Step 2. Convergence to the infinite hierarchy. In Theorem 18.11 we prove that any limit point 
Too,* = {"f£]t}k>i G 0fc>i c ([0,T],Cl) of FN,t = {lN, t }k>l witn respect to the product topology 
Tprod is a solution of the infinite hierarchy of integral equations (k = 1,2,...) 

7 <S= W( fc )(t)7So-8"«oE f ^ k \t-s)Tr k+1 \s( Xj - x k+1 ),^] (4.4) 

with initial data 7^ = (where S ^(M 3 ) has been introduced in (|3.5p ). Here Tr^i 

denotes the partial trace over the (k + l)-th particle, and U^ k \t) is the free evolution, whose action 
on /c-particle density matrices is given by 

:= e ft £i=i A . 7 ( fe ) e -«^= 1 A j (4 5) 

We remark next that the family of factorized densities, 

7 « = i^x^r, (4.6) 

is a solution of the infinite hierarchy (|4.4p if y?t is the solution of the nonlinear Gross-Pitaevskii 
equation (|3.3f) with initial data ift=o = <p. The nonlinear Schrodinger equation (|3.3f) is well-posed in 
i? 1 (M 3 ) and it conserves the energy, £(tp) := \ J \Vip\ 2 + 47rao / M 4 . From cp e fl^M 3 ), we thus 
obtain that G i? 1 ^ 3 ) for 

every i € M, with a uniformly bounded -norm. Therefore 

Tr (1 - Ai) ... (1 - A k )\ip t )(<p t \® k < \\<p t \\ k Hl < C k (4.7) 

for all t £ I, and a constant C only depending on the i? 1 -norm of ip. For the well-posedness 
of the subcritical nonlinear Schrodinger equation (|3.3|) in H , see e.g. [21]. We remark that the 
well-posedness has been established even for the critical (quintic) nonlinear Schrodinger equation in 
[151 f!6| [28] for small data and in |6] for large data. 

Step 3. Uniqueness of the solution to the infinite hierarchy. In Section 9 of [11] we proved the 
following theorem, which states the uniqueness of solution to the infinite hierarchy (14. 4p in the space 
of densities satisfying the a priori bound (|4.3p . The proof of this theorem is based on a diagrammatic 
expansion of the solution of (14. 4p . Remark that the uniqueness of the infinite hierarchy in a different 
space of densities was proven in [22J. 

Theorem 4.1. [Theorem 9.1 of [11]] Suppose T = {7 (fc) }fc>i G Q) k>1 C\ is such that 

Tr (1 - Ai) ... (1 - A fc ) 7 W < C k . (4.8) 

Then, for any fixed T > 0, there exists at most one solution T t = {7t }fc>i G ©fc>i C([0, T],C\) of 
\4-4\) such that 

Tr (1 - Ax) ... (1 - A fc ) 7t (fc) < C k (4.9) 
for all t £ [0, T] and for all k > 1. 

Step 4- Conclusion of the proof. From Step 2 and Step 3 it follows that the sequence T^t = 

{iNuk>l G ©fc>i C([0, T], C\) is convergent with respect to the product topology r pro d; in fact 
a compact sequence with only one limit point is always convergent. Since the family of densities 
T t = {7^ }fc>i defined in (|4.6p satisfies (|4.7j) and it is a solution of (|4.4p . it follows that r^t — ► Tt 
w.r.t. the topology r pro d- In particular this implies that, for every fixed k > 1, and i £ [0, T], 

7^^ — ► \(pt)((ft\® k with respect to the weak* topology of C\, and thus, by a standard argument, also 
in the trace-norm topology. This completes the proof of Theorem 13.21 □ 
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Next we prove Theorem 13.1) to this end we have to combine Theorem 13.21 with an approxi- 
mation argument for the initial A-particle wave function, which is needed to make sure that the 
energy condition (|3.4p is satisfied. This argument was already used in [p3]; we present it here for 
completeness. 



Proof of Theorem \3.1\ We assume here that, as in Theorem 13.21 the interaction potential V is such 
that |V Q F(x)| < C for all multi- indices a with \a\ < 2. We show how to remove this condition in 
Appendix [Bl 

Fix k > and x £ Cfi°(R), with < x < 1, x(s) = 1, for < s < 1, and x(s) = if s > 2. We 
define the regularized initial wave function 

~ x(kH n /N)iP n 
* N - \\x{kH n /N)^ n \Y 

and we denote by ipN,t the solution of the Schrodinger equation fjl .3j) with initial data ippf- Denote 

by T~N,t = {lNWk=i ^ e f am ily °f marginal densities associated with ipN,t- By convention, we set 

~(k) 

iNt := if A; > N. The tilde in the notation indicates the dependence on the cutoff parameter k. 
In Proposition 19.14 part i), we prove that 



k i\Tk 



(4.10) 



if k > is sufficiently small (the constant C depends on k). Moreover, in part iii) of Proposition [9TT 
we show that, for every J^> E JC k , 



WW\ m ) ^0 (4.11) 

as N — > oo. From (|4.10p and (|4.11j) . the assumptions (|3.4p and (|3.5p of Theorem 13.21 are satisfied 
by the regularized wave function ippf and by the regularized marginal densities ") Nt - Applying 
Theorem 13.21 we obtain that, for every t E R and k > 1, 



~(fc) 
Tiv,t 



(4.12) 



where ^ is the solution of (|3.3p with initial data <^t=o = 

ft) 

It remains to prove that the densities "f N t associated with the original wave function ipN,t (without 

cutoff k) converge and have the same limit as the regularized densities 7^1 ■ This follows from 
Proposition 19-H part ii), where we prove that 

\\ipN,t - i>N,t\\ = W^n - Vwll < Ck X I 2 , 

for a constant C independent of N and k. This implies that, for every we have 

, t $$\<C\\JV>\\£I*. (4.13) 
Therefore, for every fixed jfe > 1, t E R, J (fc) E K k , we have 



TrJ«( 7 $fl 



Tr 



(k) 

r N 't 



\<pt)(<pi 



< 



Trj( fc )f 7 ^ 



~(fc) 



+ 



Tr J( fc ) f^l 



l®Jfe 



< C || JWIIk 1 ^ + | T r (7$ - bt)(^r fc 



(4.14) 



Since k > was arbitrary, it follows from (|4.12p that the l.h.s. of (|4.14p converges to zero as N — > 00. 
This implies that, for arbitrary k > 1 and £ E R, 7^ — > |(/?t)(92j|® fc in the weak* topology of C\, and 
thus also in the trace- norm topology. This completes the proof of Theorem 13.21 □ 
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(k) 

5 The wave operator and a-priori bounds on j K N J t 

(k) 

In order to derive a-priori bounds for the marginal densities j N t , we need to introduce wave operators. 
We denote by W and Wn the wave operators associated with the one-particle Hamiltonian f) = 
-A+ (l/2)V(x) and, respectively, t) N = -A + (l/2)V N (x), with V N (x) = N 2 V(Nx). The existence 
of these wave operators and their most important properties are stated in the following proposition 
(we denote by s — lim the limit in the strong operator topology). 

Proposition 5.1. Suppose V > 0, with V G L 1 (R 3 ). Then: 

i) (Existence of the wave operator). The limit 

W = s - lim e i¥ e iM 

f^oo 

exists. 

ii) (Completeness of the wave operator). W is a unitary operator on L 2 (M 3 ) with 

W* = W~ x = s - lim e - iAt e~ i¥ . 

Hi) (Intertwining relations). On D{\)) = D(—A), we have 

W*t)W = -& (5-1) 



iv) (Yajima's bounds). Suppose moreover that V(x) < C(x) a , for some o > 5. Then, for every 
1 < p < oo, W and W* map LP(M. 3 ) into LP(M. 3 ), and 



\\W\\lp->lp < oo for all 1 < p < oo/, . 
v) (Rescaled wave operator). If t)N = — A + (1/2)Vn(x), with Vn{x) = N 2 V(Nx), then the limit 



W N = s - lim e ii>Nt e iAt 

exists and it defines a unitary operator Wn on L 2 (R 3 ) with 

W* N = W N X = s - lim e - iM e- if)rft . 

The wave operator Wn satisfies the intertwining relations 

W N t) N W N = -A. 

Moreover, the kernel of Wn is given by 

W N (x; y) = N 3 W{Nx; Ny) and W N (x; y) = N 3 W*(Nx; Ny) 

where W(x;y) and W*(x;y) denote the kernels ofW and W* . In particular, it follows that, if 
for every 1 < p < oo, the norms 

\\Wn\\lp^lp = \\W\\lp-^lp < oo and \\W n \\lp->lp = \\W*\\lp^lp < 00 
are finite and independent of N. 
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Proof. The proof of i), ii), and iii) can be found in [24]. Part iv) is proven in [30|, 129]. Part v) follows 
by simple scaling arguments. □ 

In the following we will denote by Wnj\ and, respectively, by Wj^Uj\, the wave operators W and 
Wm acting only on the relative variable Xj — X{. In other words, the action of Wuj\ on a iV-particle 
wave function ifjjy S L 2 (M. 3N ) is given by 

dv W(xj - xf, v) t/) N I xi,..., 1 2 3 +-,..., * 2 x N \ (5.2) 

if j < i (the formula for i > j is similar). Here W(x;y) is the kernel of the wave operator W. An 
analogous formula holds for the rescaled wave operator Wn- Similarly, we define W% and ,.^y 
Using the wave operators, we have the following energy estimate. 

Proposition 5.2. Suppose V > 0, V G L 1 (R 3 ) and V{x) = V(— x) for all x 6 K 3 . Then we have, 
for every i ^ j, 

(^ N ,H%^ N ) > CN 2 J dx |(Vi • V,-) 2 W^ {t ^ N \ 2 , (5.3) 

where ^ ^ denotes the wave operator defined in Proposition I5.il acting on the variable v = 
Xj — Xi (defined similarly to 15. Sty ). 



Proof. We define, for j = 1, . . . , N, 



Aj + ^^VNixi-Xj). 



2 . 



Then we have Hn = Ylf=i hj an d thus 
(i> N ,H 2 N i> N ) >N{N-l)(i> N Mh2i>N) 

= N(N - 1) ( ^v, -A x + £ V V N ( Xi - xi) ) I -A 2 + ~ V V N ( Xj - x 2 ) I 



^-Ai + i^Vjv(xi-a:i)j ^-A 2 + i^VW( a 



> iV(iV - 1) ( Vjv, ( -Ai + ^^v(x! - a? 2 ) ) I -A 2 + ^(a* - x 2 ) ) rj) N 



(5.4) 



Now we define the new variables 

21 + 22 



and w = xi — 22 • 



Then we have 



and thus 



We set 



Vi = -V« + V„ and V 2 = -V u - V t 



-A u + A v + V u - V„, and A 2 = yA„ + A„ - V„ • V„ 



h v = -A v + -Vn(v) 
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Then 



■ IInVn) > N(N-l)(i> N , ( --A u + h v + V u - V v j (--A u + h v -V u -V v ) r v 



N{N - 1) ( 



1 \ 2 1 

-jA u + h v \ -(V u -Vv) 2 +2 V u-(VV N (v)) 



(5.5) 

IpN ) ■ 



Next we note that 

(5.6) 

dudv ij)ff(u + v/2,u — v/2,Xn-2)Wn(v) • V u iPn(u + v/2, u — w/2,xat_ 2 ) = . 

In fact, by the permutation symmetry, iJ)n(%i, x 2, x iV-2) = (^2, xjv-2)- This implies, in the 
it, f-coordinates, that tp]y(u+v/2, u—v/2, Xj\r_ 2 ) = '4 , n(u—v/2, u+v/2, xjv-2) an d also that V u ^n(u+ 
v/2, u — v/2, X7V-2) = ^u^n{u — v/2, u + v/2, xjv_2). On the other hand (VV/v) (— i>) = — (VVjy) (f ). 
Therefore, the integrand in (|5.6p is antisymmetric w.r.t. the change of variables u — ► —v, and the 
integral vanishes. 

Using also that 

(V„ • V,,) 2 < (-A„) (-A^) < (-A u ) ^ , 

it follows from (15.51) that 



(ip N ,H 2 N ip N ) >N(N-l)(ip N , [~A u -h}j ^ . (5.7) 

Next we make use of the wave operator Wjv defined in Proposition 15. 1\ acting on the variable 
v = X2 — x\. By the intertwining relations (|5.ip . we find 



(1>n, H% $ N ) > N(N - 1) (w^ {h2) ip N , Qa u - A,) H^ j(1>2) ^\ . 



(5.8) 



In terms of the coordinates x\ and X2, we have Vi-V 2 = (1/4)A U — A„. Therefore, by the permutation 
symmetry, the bound (|5.8p implies (|5.3p . □ 

Proposition 15.21 implies strong a-priori bounds on the solution of the iV-particle Schrodinger 
equation. 

Proposition 5.3. Suppose that V > 0, V £ L 1 (M 3 ), and V(—x) = V(x) for all x G M 3 . Let tfj N j be 
the solution of the Schrodinger equation 111.3]) , with initial data satisfying the assumption fl^.^l ) (with 
k = 2) of Theorem \3.3l and let {j^\}^ =l be the marginals associated with ipN,t- Then, for every 
1 < j < N, we have 

(i/>N,t , (1 - &j)i>N,t) < C and thus Tr(l- Aj ) 7 jg < C 

for every 1 < j < k < N ( and for a constant C which only depends on the initial data ipN through 
the constant on the r.h.s. of [3.$ )- Moreover, for any i ^ j, 
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uniformly in N > 1 and in t £ M. i/ere Wjv,(ij) denotes the wave operator Wn defined in Proposi- 
tion \5.1\ acting on the variable Xj — X{. In terms of density matrices, we obtain the a-priori bounds 

Tr ((V, • V,) 2 - A, - Aj + 1) W* m>j) 4} t W N>[i>j) < C 

uniformly in N > 1 and in t £ R and /or all 1 < i < j < k (with a slight abuse of notation, we 
denote here by W^^fi and ^ the operators acting on the k-particle space L 2 (M? k )). 

Proof. The first bound follows simply by the symmetry of the wave function, by energy conservation, 
and by the condition V > 0. To prove the second bound, we compute 

W,^, ((V, • V,) 2 - Aj - Aj + l)W* Nm ^ N , t ) 

= (WN,(i^N,t, (Vi • V 2 ) 2 W^ (lj2) Viv,t) (5.9) 
+ W,(i )2 )^,t, (-Ai - A 2 ) W* Nt{l ^ Njt ) + 1 • 
The first term on the r.h.s. of the last equation can be bounded by 

W,(i,2)^,t, (Vi • V 2 ) 2 W^ )(1)2) ^,t> < CN- 2 (if> Nt t,Hfrl, N ,t) 

= CN~ 2 (ip Nfi ,H 2 N ^N,o) < C 
using Proposition 15.21 and ([3.4|) . The second term on the r.h.s. of (|5.9p is estimated by 

( W N,(l,2)^N,t, (- Ai - A 2 ) ^,(1,2)^,*) 

= 2 /w^ (li2) V>iv,t, f-Ao,!-^ - ^A (a . 1+a . 2 )/ 2 J W^^ 1>2 ^N,t) 

= 2 Z^JV.t, ^-A X1 _ X2 + ~Vn(xi -x 2 )-j A {xi+X2)/ ^j ip Njt \ ( 5 - n ) 

= {ipN,t, (-Ax! - A X2 + Viv(xi - x 2 )) ipN,t) 

2 2 

< j^{tpN,t,HNlpN,t) = J-f(lpNfl,H N lp N fi) < C . 

□ 



(5.10) 



6 Compactness 

In this section we prove the compactness of the sequence Tnj = {7jvi}fc>i w - r -t- the topology r pro d 
(defined in Section 2]). 

Theorem 6.1. Let the assumptions of Theorem \3.S\ be satisfied and fix an arbitrary T > 0. Then 
the sequence Tnj £ ©fc>i C([0> T], C\) is compact with respect to the product topology T pro d generated 

by the metrics rj^ (defined in Section^. For any limit point T^t = {7oo\}fc>l? looi *s symmetric 
(k) 

w.r.t. permutations, 7^ > 0, and 

Tr^] t <l (6.1) 

for every k > 1. 

(k) 

Proof. By a standard argument it is enough to prove the compactness of 7^ t for fixed k > 1 with 

respect to the metric r/fc. To this end, it is enough to show the equicontinuity of 7)^ with respect 
to the metric A useful criterium for equicontinuity is given by the following lemma, whose proof 
can be found in [11, Proposition 9.2]. 
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Lemma 6.2. Fix k G N and T > 0. A sequence 7 jyl G iV = /c, fc + 1, . . with 7]^! > and 
Tr^\ = 1 /or a// i G [0,T] and N >k, is equicontinuous in C([0, T],£^,) w;ii/j respect to the metric 
i/ and on/y i/ i/iere exists a dense subset of K,^ such that for any G and for every 
e > i/iere exists a 5 > suc/i t/iai 



.,(*) 



sup 

Af>l 



< e 



(6.2) 



/or a// 1, s G [0, T] wii/i |t — ,s| < S. 

We prove (1631) for all jW G /C fc such that |||j( fc )||| < 00, where we introduced the norm 

|||JW|||= sup^) 4 ...^) 4 ^) 4 ...^) 4 



x [ |j( fc )(x' fe ;x fc )|+^(|V, j j( fc )(x' fc ;x fe )| + |V 1 , jW(x' fc ;x fc ) 
3=1 



(6.3) 



It is simple to check that the set of g for which |||J^||| < 00 is dense in /C^. 

Rewriting the BBGKY hierarchy (|2.ip in integral form and multiplying with an arbitrary observ- 
able j( fc ) G /Cfc with I I < 00, we obtain that, for any r < t, 



^ (7$ - 7^) | < E / H 1 * J{k) [- A * 7$] I + E f H 1 * J{k) \y»& - 

j=l r i<j r 



(k) 
T Ns 



3=1 



N 



E [ ds Tr J( fc ) \nV n ( Xj -x fc+1 ) j7 <$ 



(fc+i) 



(6.4) 



To control the first term on the r.h.s. of the last equation, we observe that, using the notation 
Sj = (1- Aj) 1 / 2 , we have 



Tr jW[-A il7 $] = Tr J^S, - ^ J^ 1 ) ^7$^ 



,(*) 



< 
< C 



S^J^Sj + SjJWSJ 1 jTrCl-A^S 



,(*) 



(6.5) 



Here we used that, by Proposition 15.3 



sup Tr(l - A.-hS < C 
uniformly in N. 

To bound the second term on the r.h.s. of (16. 4j) . we decompose 7 jy^ a = X^A^ |Q )(Q | for 



ef } g l 2 



with ||£^|| = 1) A^ > 0, and £^A^ = 1 (here we omitted the dependence of £^ 



and X\ on N,s from the notation). Then we find, for example for the term with i = 1, j = 2, 

TrJ^V N (x 1 -x 2 ) 1 ^ s = J24 k) J dx * dx fe J (k) «^k)V N (x l -x 2 )^ (x,)?f (x' fc ). (6 . 6) 
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Denoting by Wjv the wave operator associated with the Hamiltonian f) 



W, 



N,(i,j) 



N = -A + (l/2)V N (x), and by 

the wave operator Wn acting on the variable Xj — Xi (as defined in (|5.2p ). we can estimate 
(introducing the new variables u = (x\ + x 2 )/2 and v = x\ — X2) 

f dx k dx^ k \x' k ,x k )V N (x 1 -x 2 )if\x k )f k \x^) 

= j dudvdx 3 . . . dx k dx k J^ k \x' k ;u + v/2,u — v/2, x 3 , . . . , x k ) 

x V N {v)if ] (u + v/2, u - v/2, x 3 ,...,x k ) lf ] (4) 
< J dudvdxs . . . dx k dx' k J^(x' k ; u + v/2, u — v/2, x 3 , . . . , x k ) — J^ k \x' k ; u, u, x 3 , . . . , x k ) 
x V N (v)\^ k) (u + v/2,u- v/2,x 3 , . . .,x k )\ \Z (k \x' k )\ 
+ J dudvdx 3 ...dx k dx' k \J^ k \x' k ;u,u,x 3 ,...,x k )\\(W^V N )(v)\ 
x (w* N ^ 2 / k) )(u + v/2, U -v/2,x 3 ,...,x k )\ 

where in the last line we used the L 2 -unitarity of the wave operator in the ^-variable (before taking 
the absolute value inside the integral). Hence 



,x 3 , ..,x k )\\v\V N (v)\£ 



dx k dx' k j( k \x> k ,x k )V N (x 1 -x 2 )&\x k )& ) (x' h 
-YlJ dud-i;dx 3 • • • dx k dx' k J dr|V Xj J (fc) (x' fc ; u + 

2 f f 1 TV TV 

+ dudvdx 3 ■ ■ ■ dx fc dx' fe J dT\V Xj J^ k \x k -u + —,u- —,x 3 ,..,x k )\\v\V N {v) 

x \tf\u + v/2,u-v/2,x 3 ,...,x k )\ 2 
+ J dudvdx 3 . ..dx k dx' k \J { - k \x' k ;u,u,x 3 , . . . ,x k )\ \{W%V N ){v)\ \x' k )\ 2 

+ J dudvdx 3 ...dx k dx' k \J^(x' k ;u,u,x 3 ,...,x k )\\(W^V N )(v)\ 

x (j^,(i,2)Zej ( u + V A u - V / 2 ' x 3 ,...,x k ) 
Using the norm defined in (|6,3p . we find 

J dx k dx' k J (h \^k)V N (x 1 - x 2 )^ k \x k )f k \x' k ) 



<C fc |||j( fc )|||||e. 



(fe)l|2 



dv \v\ Vn(v) 



+ C k I I / dt;dudx 3 • • . dx N \v\V N {v) \6 k \u + v/2, u - v/2, x 3 ,...,x k )\ 



+ C k \\J^\\U i i k) \\ 2 J dv\(W* N V N )(v)\ 

+ C k \\J^\\ j dudvdx 3 ...dx k \(W* N V N )(v)\ \(w^ {li2) ^(u + v/2,u-v/2,x 3 ,...,x k ) 



(6.7) 
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Since, by scaling || \v\V N ||i < CN~ 2 and || \v\V N || 3/2 < CN" 1 and since ||W£V/v||i < C\\V N \\i < 
CN" 1 (by the Yajima's bounds in part v) of Proposition 15. 1|) . we conclude that 



j dx fc dx / fe j( fc )(x' fc ;x fc )^(x 1 -x 2 )ef(x fc )ef ) (x / A 

^ c k \jv>i 



N 



W* NA i,2)&\ ((Vi • V 2 ) 2 - A x - A 2 + 1) W* N){1)2) t 



(6.8) 



for a /c-dependent constant C k . Here we used Lemma 1 10 . 1 1 to bound the last term on the r.h.s. of 
(|6.7p . and that, by the Sobolev inequality H^^H^e < C , ||V„^ fc ^||^2, 

dudx 3 . . .dx k dv \v\V N (v)\^ e k \u + v/2,u - v/2,x 3 , . . . ,x k )\ 2 

< \\\v\V N \\ 3 / 2 J dudx 3 ...dx k dv\V v ^ k \u + 



<CN- l Uf\(-^)i 



(k) 



(6.9) 



< CN' 1 (tf\ (-A, + (1/2) VM) ef) 

< CAT 1 (w-* )(1)2) ^ (fc) , (-Aj - A 2 ) ^ )(1)2) # ) 

to bound the second term on the r.h.s. of (|6.7p (applying the intertwining relations (|5.ip and adding 
a positive term — A n as in (15. 111) ). From (16.61) . (16. 8p . and from Proposition 15.31 we obtain that 



Tr J( fe ) V)v(xi - x 2 ) 7 



(AO 
JV,s 



< Tr ((Va • V 2 ) 2 - A : - A 2 + l) ^tJXm 



< 



N 

g fc |||^ (fc) lll 
iV 



(6.10) 



for all s G R and a constant Cfc only depending on k (and on the constant appearing on the r.h.s. of 
(|3.4p ). Similarly to (|6.10p . we can also show that 



Trj( k hP 9 V N (x 1 -x 2 ) 



< 



C k \\J^\\ 



N 



(6.11) 



Since (|6.10p and (|6.1ip remain valid for all summands in the second term on the r.h.s. of (|6.4p . we 
obtain that, for all fceN, for all t G [0,T] and for all j( fc ) G fC k with |||j( fc )||| < oo, 



J^ \Tr [Vtfixi-xj),^ 

i<j 



< 



Ci- 



N 



(6.12) 



for all seR. 



Also the third term on the r.h.s. of (|6.4p can be bounded similarly. In fact, using again the 
decomposition = IC^* )(£i fc+1 ^l we have, for example considering the term with 

J 1- 



TvJ^NVNixx-Xk^)^ 



A? +1) 1 dx fc dx' fc dx fc+1 J«(x / fc ;x fc )iV^(x 1 - x fc+ i)£? +1) (x fc , x k+1 )& +1) (x' fc , 



Xfc+i) • 



(6.13) 
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The absolute value of the £-th summand can be estimated by 

J dx fe dx'^dxfc+i jW(x' fe ;x fc )iVViv(a;i - x k+1 )if +1) (x fc , x k+1 )lf +1 \x k , x k+1 ) 
< J dudvdx 2 ... dx k dx k NVn(v) \^ k+1 \u + v/2, x 2 , ... ,x k ,u — v/2)\ 
x j( fc ) (4; u + v/2, x 2 , . . . , x k )lf +1) (4 , u - v/2) 

- j(k) (4 ; u > x 2 , ■ ■ ■ , x k (4 1 u ) 

+ y dMdudx 2 ...dx Jfc dx / fc iV|(W^Viv)(w)| \J {k) (x' k ;u,x 2 , . . . ,x k )\ 

x |ef +1) (x' fe , u)| | (w* NXl>k+1) & +1) ) (u + v/2, x 2 ,...,x k ,u- v/2) 



I + II. 



Here 



dudi;dx 2 • • • dx k d~x' k NVn(v) \t;i k+1 \u + v/2, x 2 , ■ ■ ■ , x^, u — v/2)\ 



' d 
dT d7 



< / diidudx2 • • • dx k d~x.' k I dr NVn(v)\v\ \Q k {u + u/2, x 2 , ■ ■ ■ , x k , u — v/2)\ 



"jfc> u 7T , 







ru, 



+ 1 ^ (fc) (4; u + y > x 2 , ■ ■ ■ , x k ) 1 1 v Xk+1 $ +1) (4 , « - ~2 



ru , 



TV , 



(6.14) 



(6.15) 



Through a weighted Schwarz inequality, we find 

I < J dudvdx 2 . . . dxfedx' fe J dr NV N (v)\v\ \V Xl J (fc) (x' fc ; u + x 2 , . . . , x k )\ 

x (iV-V2 + w/2> x 2 , . . . , x fc , u - ^/2)| 2 + ArV2|^+i) ( 4 ; 

+ y dudf dx 2 . . . dx k dx' k J dr NVn(v)\v \ \ (x' fc ; u + x 2 , . . . , x k )\ 

x |^i) (t| + w/2| x 2 j ••• j Xfc, u <;/2)| 2 + iV 1 ^ I Vxfc+I ef +1) (4, „ - ^)| 2 ) . 

Extracting the observable from the integral (after integrating some of its variables and taking the 
supremum over the other), and using Sobolev inequalities where needed, we find 



I < C k iV 1/2 ||MF7v||3/2 II J(k) \l J dudvdx 2 . . . dx k \V v i ik+l \u + v/2,x 2 ,...,x k ,u-v/2)\ 2 

+ C fc iV 3 / 2 |||,;|VW|| 1 |||j( fc )||| || 2 + / dudx' k \V k+ 4 k+1 \x' k ,u)\A (6.16) 



< 



N l/2 \ W N,(l,k+l)Q » l- A l - A ^+l + ^(l.fe+l)^ /" 



In the last line we proceeded similarly to (|6.9|) for the two terms with derivatives, using the bounds 
IIMVarlls/a < CN- 1 and |||«|Vat||i < CN~ 2 . 
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As for the second term on the r.h.s. of (|6.14p . we can bound it applying a Schwarz inequality 
and Lemma ll 0.1 1 by 

II < J dudvdx 2 . . .dx k dx! k: N\{W* N V N ){v)\ \J^(x' k ;u,x 2 , • • • , x k )\\$ +1) (x' fc , u)\ 2 
+ j dudvdx2...dx k dx' k N\(W^V N )(v)\\J {k) (x' k -,u,X2,...,x k )\ 

x ( W N,(i,k+i)& +1) ) (u + v/2,x 2 ,...,x k ,u-v/2)\ 2 

< c k N\\j^\\ \\wzv N \\i (w* N ^ k+1) ^\ [(V! • v k+l f - Ai - A k+1 + 1] w* Ni(1;k+1) ^ k } ) . 

Since, by Yajima's bounds (Proposition 15.11 part v)), ||W^yAr||i < C||Vw||i < CN~ l , it follows that 

II < C k I J( fc )||| (w* N ^ k+1) if\ [(Vi • V fe+1 ) 2 - Ax - A fc+1 + 1] W* Nt{hk+1) & } ) . 
Inserting this and (|6.16p into the r.h.s. of (|6.14p . it follows from ()6.13j) after resumming over I that 
TrJ^NVNixt-Xk+ih^ 



u ( fc + 1 ) 



< C k I JW I Tr ((V! • V fc+1 ) 2 - Ai - A fc+1 + l) W^^X'W^ k+s } 
<C k 



for all s G K (in the last line we used Proposition 15. Since the same bound remain valid if we 
replace x\ with an arbitrary Xj, j = 2, . . . k (and also if the potential lies on the right of the marginal 
density 7^^), it follows that 



£ |tv j« [nv n ( Xj - x k+1 ),^ 1] ] | < C k I J« I . 
i=i 



(6.17) 



From dS3D, dSISD, (16TT21) . and (I6TT7]) . it follows that 

TrJ« {lP >t -lf^\<C k \\jM\\\t 



s . 



This implies (|6.2p . and thus the equicontinuity of the sequence T^t = {7yvl}fcLi with respect to the 
metric r proc j. 

The proof of the fact that 7^ t is symmetric w.r.t. permutations, that it is non-negative and such 
that Tr7^ t < 1 can be found in [I2~l Theorem 6.1]. □ 



7 Higher order a-priori estimates on the limit points T^j. 

The goal of this section is to establish strong a-priori bounds for the limit points Toot of the sequence 



N,t 



{7tv i}fc=i • ^ s we m [12]' we wm obtain strong a-priori estimates on r M| ( by proving higher 
order energy estimates, which compare the expectation of powers of the Hamiltonian (ip^, H^i/in) 
with certain Sobolev norms of the iV-particle wave function ipN . It turns out that the expectation of 
powers of the Hamiltonian can only control the Sobolev norms of ipN i n appropriate regions of the 
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configuration space. To characterize these regions, we introduce the same cutoffs we used in |12| . 
For a given iength scale £ > (in our analysis, we will need that iV -1 / 2 <C £ <C A^ 1 / 3 ), we set 



h(x):=e 1 — . (7.1) 

Note that h — if \x\ ^> £, and h ~ e _1 if |x| <C £. For i = 1, . . . , N we define the cutoff function 

0i(x) := exp ^-1 ^ h( Xi - Xj )\ (7.2) 



for some e > 0. Note that 0j(x) is exponentially small if there is at least one other particle at distance 
of order £ from Xj, while #i(x) is exponentially close to 1 if there is no other particle near xi (on the 
length scale £). Next we define 

0< n >(x) := ^(x) 2 " = exp ^ £ h ^ ~ X ^J t 7 ' 3 ) 
and their cumulative versions, for n, k G N, 

ejW(x) := 0? } (x) . . . *W(x) = exp ( £ £ fcfc - x,) ) . (7.4) 

(n) 

To cover all cases in one formula, we introduce the notation Q k = 1 for any k < 0, n G N. Some 

important properties of the function 0jj. , used throughout the proof of Proposition ^. 11 are collected, 
for completeness, in Lemma lA.ll 

Proposition 7.1. Suppose that V > 0, mi/i |V Q y(x)| < C /or a// |a| < 2. Let V G L 2 (M 3Ar ) be a 
function symmetric in all its variables. Suppose that £ S> N~ l l 2 (in the sense that there exists 5 > 
with N 1 / 2 £ > N s ). There exists Cq > such that for every integer k > 1 i/zere exists Nq = No(k) 
such that 

^,(H N + Nfy)>C%N k J dx6^ 1 (x)|V 1 ...V^(x)| 2 (7.5) 

for all N>N . 

Proof. We use induction over k. For k = 1 the statement follows directly from Vjv > 0, since on the 
symmetric subspace 

TV 

H N + N > V*Vi = ATV*^- (7.6) 
i=i 

for any fixed j = 1, 2, . . . , N. We will present the k = 2 case in details and then comment on the 
general case. Set T = Hn + N > for brevity and use the induction hypothesis 

T 2 > C A^T 1 / 2 V^ViT 1 / 2 

> CoNT^^iefVtT 1 / 2 

> ^C Q NVlT l l 2 e\T 1 / 2 V l - CoNiT 1 / 2 , Vifflfrr 1 / 2 , Vi] ^ 

> -c Arv^?re 2 Vi - CoAfv^r 1 / 2 ,^ 2 ]*^ 1 / 2 ,^^! - c n[t^ 2 , Vip^T 1 / 2 , Vi]. 
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In the first term we use that Hn + N > Yl!j=2 to obtain 

-CoNVlefTOlVi > ]c N(N - l)V^ 2 V^V 2 2 Vi 

4 \ ' (7-8) 

> -C iV 2 V^V^V2Vi-CoiV 2 vnV2^i 2 ]*[V2^i 2 ]Vi 

for all N large enough. Since 9^ = 9f, we would obtain ([73]) for k = 2 with C < 1/8 once we 
show that the commutator terms in (|7.7p and (17. 8h are negligible. 

The commutator in (|7.8p on symmetric functions can be estimated by 



r n 2 

C ^ 2 Vj[V 2) e?]*[V 2 ,5?]Vi = ^-^VKV^ 2 )^! 

i=2 

< o(r 2 iv- 1 )iv 2 ViVi < o(r 2 iY _1 )r 2 = o(i)r 2 

where we used (jA.ip . recalling that 0\ = 0^ , and we also used T > N and (|7.6p . 
To estimate the two commutators in (j7.7h . we express 



(7.9) 



1 f 00 1 1 
[T 1 / 2 , A] = - / — - [A, T) — - s l ' 2 ds (7.10) 
7T J T + s T + s 



(7.11) 



for any operator A. 

To estimate the first commutator term in (j7.7|) we note that, by Schwarz inequality, 
[r i/ 2) ^ [T i/ 2) ^ < c (lQg K) ^ ^ rr _L_ [02) T] (s)2 ds 

+ C / - [9l T}* — ^ [fl 2 , T] (s) 5/2 ds, 

7 K T + s 1 1 J (T + s) 2 L x ' J T + s w 

where if = exp(iV e ) for some e > 0. Estimating (T + s)~ 2 < (s)~ 2 (using T > iV), we have 

r x i i 

iWi[T 1 ' 2 ,0 2 ]*[T 1 ' 2 ,0 2 ]Vi < cN 1+E / VI- [T,# 2 ]*[T,# 2 ] V x ds 

Jo r + s r + s , . 

/■OO 1 1 v ' 

+ civ / vi-—[T,e 2 1 r[T,e 2 l ]-—v 1 (s) 1 / 2 d s , 

J K T + s T + s 

and we can estimate 
[T^efme 2 ] = £ (2V* • (V^ 2 ) + (A^ 2 ))(2(V^ 2 ) • V, + (A^ 2 )) 

< c]T [v* • (V^ 2 )(V^ 2 ) • Vi + V* • |A^ 2 |V, + |V^ 2 ||A^ 2 ||V^ 2 | + |A^ 2 ||A^ 2 

(7.13) 



22 



by using (|A.ip and (|A.2p and the fact that 9 2 = 0j . Thus, the first commutator term in (|7.7p is 
estimated as 

ivvj [r 1 / 2 , e 2 }* [t 1 / 2 , ^ 2 ]Vi < o(iv 1+e r 2 ) / vj — !— r — !— Vids 

./() T + s T + s 

/■OO 1 1 

+ 0(N£- 2 ) / VJ=— T— - ViM^d* 

< o(iv 1+e r 2 ) v*Vi + o(Nr 2 K~ 1/2 )v* t Vi 

< o(iv £ r 2 )r + o(i)r 2 

< o(l)T 2 

if we choose e > so small that l~ 2 <C iV 1-5 . When estimating the term V^TVi in the last step, 
we could afford estimating any commutators, since K~ 1 / 2 is exponentially small: 

v ^ TV i = jf E V ,* TV ^ = ^ E [(- A ^ T - v ,* ■ E( vv ^)(^ - * 



iV^T 2 - AT' 



- 1 ^{xi - Xj)T - V* • (Vy 7 v)(x i - x,) 



< 2N-'T Z + 0{N*) 

using that |Vjv(x)| < CN 2 , and \W N (x)\ < CN 3 , for all x £ M 3 . 

Finally, we estimate the second commutator term in (|7.7p by using again Schwarz inequality in 
(I7.10P but this time we do not split the integration: 

[T^ 2 ,Vi]*&t[T^ 2 ,Vi] 

l-CO 1 11 1 

Jo T + s T + s T + s T + s (7.15) 

1-tX) 1 11 1 

^ ciV E / ^^V N )(x l -x l )-^9f-—(VV N )(x 1 -x l )-—(sf/ 2 d S 
f—' y T + s T + s T + s T + s 



where we used [Vi,T] = X^i(VV/v)(xi ~~ 

Since T = £^ ■ — Aj + iV is a positivity preserving operator and V > 0, T is also positivity 
preserving and its resolvent kernel satisfies 

1 / 1 / 

(x;y) < - (x;y), 



T + s 



Tn + s 



and thus 



I : = 



(VVtv)(xi 



' * 4 ' -(VVat)(xi — Xj 



T+s "T+s 



< 



KWjvXxi-xi)! 



i 



T + s 



-4£- £ h(x 1 -x i ) 



Tn + S 



|(VFiv)(xi-Xi)| 



(7.16) 



where we also estimated 9\ by keeping only one summand in its definition (|7.2p . Introducing the 
variable y = x\—Xi, and observing that T 2 (R 37V , dx) ~ L 2 (R 3 , dy; L^IR 3 ^ -1 ), dzdx2 . . . dxj . . . dxjv)) 
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(where the hat means that the variable Xi is omitted), we obtain that 

1 —Ae-ehf„A 1 



I < sup 

M>0 

< sup 

M>0 



l(VVW)(y)| 
\(yV N )(y)\- 



-A v + M + N + s 



-4£-^h(y)_ 



A v + M + N + s 



\(yV N )(y)\ 



1 



+ M + N + s 



-ii- s h{y)_ 



1 



A y + M + N + s 



HS 



where the norms on the last two lines are, respectively, the operator norm and the Hilbert-Schmidt 
norm of an operator over L 2 (R 3 ,dy). The last equation implies that 



/ < / dydy'e- 41 "^) 
/ -4£- e ft(») 



1 



'M2 



|VVKr(i/)| 



p -2 v / 7V+i|?/- 3 /'| 

< / dydy'e-« ^ - | V Viv(y' 



/\|2 



(7.17) 



< 0(e" 



since /i ~ e 1 on the support of Vat. We will use this bound for s < K := exp(c£ £ ) with a sufficiently 
small c > 0. From (|7.16p . we also have the trivial bound 



I < 



N 6 
W 



(7.18) 



that will be used for large s. Inserting these estimates into (|7.15p . we have 



r K 

[T 1 / 2 , Vi]*^^ [T 1 / 2 , Vi] <0{N 2 e' rE ) / 

Jo 



(s) 5 / 2 ds 
(T + sf 



0(N 8 



K 



{T + sf 



0(e 



-cl~ 



i.e. this commutator term is sub exponentially small in N and this completes the proof of (|7.5p for 
k = 2. 

The proof for general k > 2 follows the same pattern as for k = 2. Introduce the notation 

D k = ViV 2 ...V fc . 
We recall the summation convention: for any operator A, we denote 

3 3 

DlAD k := V . . . V V* . . . V* AV Xt . . . V X1 

Ol = l Ofe = l 

where Xj = {xj t \,Xj : 2,Xj$) are the three coordinates of Xj € R 3 . 

Using the induction hypothesis, Bj._ 1 > [©j. ] 2 = ®fc and (|T.6j) we obtain, similarly to (|7.7|) 
and (USD, 



(7.19) 



T k+1 > C^T^Dlefl^T 1 ' 2 

> C^N k T^ 2 Dl[e { k k) ] 2 D k T^ 2 

> \c k ,N^Dl +l @t +1) D k+l -C k Q N k Dl[^ 

_ c fc 7V fe [r 1 / 2 , J D fc ]*ef +1) [r 1 / 2 jJ D fc ] _ c k N k+1 D* k [v k+1 ,e {k) }*[v k+1 ,e k k) ]D k , 
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for all N sufficiently large (depending on k). The first term gives the desired result if Co < 1/8; in 
the sequel we show that all three commutator terms are negligible. 

The first commutator in (|7.19|) is estimated exactly as the first commutator in (|7.7p . after replacing 
Q\ = @i with s[ ■ The estimates (|7.11[) are ()7.12j) are identical for k > 1 as well. In the key estimate 

(I7.13p . the only properties we used about 9f = 0^ from Lemma fA.ll were those that hold for oj^ 
as well. 

The last commutator in (17.190 can be estimated similarly to (|7.9p by using (lA.lj) 

Ik ATk+l N 



j=k+i 

< o{r 2 N- l )N k+l D* k ® { l:~ l) D k 

< O k (r 2 N- 1 )NT k = o k {l)T k+1 



(7.20) 



by the induction hypothesis and T > N (here we use the notation / = o k (g) if f j g —* as N — > oo 
for fixed k; analogously for / = O k (g)). 

Finally, the estimate of the second commutator in (|7.19p is similar to that of the second com- 
mutator in (|7.7j> . but more commutators need to be computed. Similarly to (17.150 and taking the 
permutation symmetry into account, we have 

[T^,D h YQ^ +1 \T l /\D k \ 

poo -| 1/1 1 

< C kNj2 / Tfr— DU(VV N )(x k -x i ) 7 ^e[ k+1) 7 ^(VV N )(x k -x i )D k . l7 ^(s) 5 / 2 ds 
Jo J + s i -\- s 1 -j- S 1 -h s 

+ c k / -—D^ 2 (v 2 v N )( Xk - Xk ^)-—e k k+1) 

Jo 1 + s 1 i s 

x 7 ^—{V 2 V N ){x k - x fc _ 1 ) J D fe _ 2 -^( s ) 5 / 2 ds 
T + s T + s 

(7.21) 

We will need the following lemma whose proof is postponed. 

Lemma 7.2. Let ip £ L^(IR 3Ar ) be a function symmetric in all its variables and let 5 > 0. Choose a 
strictly increasing sequence of positive constants {c k } k >i. Then for every integer k > 1 there exists 
Nq = No(k,S) such that 

{H N + N) k ip) > e~ CkN& J dx |Vi . . . V fe V(x)| 2 (7.22) 

for all N > N . 

We demonstrate the estimate of the first term in ()7.2ip . the second one is similar. Using & k k+1 ^ < 
e -i~ e h(x k -xi) ) we b tain) s i m il ar ly to ([7TTH|i - ([7TTH]> that 



(VV N )(x k - Xl )^-Q k k+1 ^^—(VV N )(x k - Xl 



< 0(e- e ~ ) 



and also 

w 
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Let K := exp(c£ £ ) with a sufficiently small c > 0. Choosing a sufficiently small 5, so that N s <C £ £ , 
by using (|7.22p . we have, 

Jo [T + s) 2 J K {T + sy(s} 2 

< 0{e- c '^ e )T k - 1 < o{l)T k+1 . 

This completes the proof of Proposition 17.11 □ 



Proof of Lemma \ 7. %\ We proceed by a step two induction on k; for k = 1 the claim follows from 
(|7.6p . We now consider the k = 2 case. Similarly to (|5.4p . but keeping also the /i 2 terms in the 
expansion of -ff 2 ^, we find 

T 2 > N(N - 1)( - Ai + ^jv(xi - s 2 )) ( - A 2 + ivjv(xi - x 2 )) + J\r( - A x + - ^ 2 )) 2 

> (JV 2 /2)(l>p2 - 2V^V! - 2V^V 2 - 4|| WtvIIL) + ^((V^) 2 - 2V* 1 V 1 - 4||V^||L) 

> (N 2 /2)D* 2 D 2 + iV(ViVi) 2 - CiV 2 V^Vi - CN S 

> (N 2 /2)D* 2 D 2 - 0(N S ). 

(7.24) 

Combining this bound with T 2 > N 2 , it follows that 

T 2 > cN- A D* 2 D 2 (7.25) 

for a sufficiently small positive c. 

Now we show how to go from k to k + 2. By the induction hypothesis, we have 

T k+2 > e -c k N s TD * DkT > e -c fc iV* (±. D * T 2 Dk _ 2 [£) fc , T]* [Z) fc , T]) . (7.26) 

In the first term we can use (1735]) in the form T 2 > cN~ i V* k+1 Vl +2 V k+2 V k+1 , which holds for all N 
large enough (because of the factors D k , we only have symmetry on the last N—k variables; this means 
that instead of (17341 . we are going to obtain T 2 > (TV- k)(N-k- l)A k+1 A k+2 > (N 2 /2)A k+1 A k+2 
for all N large enough). 

The commutator term, after several Schwarz inequalities, can be estimated as 



[D k ,T]*[D k ,T] < C fc (iV 2 ^,_ 1 ||Vy 7 v|| 2 J Dfc-i + ^.-2l|V 2 ^v|| 2 ^- 

< C k N 8 (l^l + D* k _ 2 D k _ 2 ) (7.27) 

< C k N 8 e Ck -^ NS T k - 1 < C k N 8 e c ^ NS T k+2 , 

where we used the induction hypothesis for k — 1 and k — 2 and, by convention, D m = 1 for m < 0. 
Inserting this estimate into (|7.26p . we obtain 

T k+2 > cN -A e - Ck N s D * k+2 D k+2 - Ck N 8 e-^- Ck -^ NS T k+2 . 

Since c k is strictly increasing, we obtain (|7.22p for k + 2. 

Actually the proof shows that a sufficiently large ^-dependent negative power, N~@ k , would suffice 
on the r.h.s. of (|7.22[) instead of the sub exponentially small prefactor. □ 
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The higher order energy estimates proved in Proposition 17. II are used to show the following strong 
a-priori estimates on the limit points = {7^ t }k>i of the sequence r^t- 

Theorem 7.3. Suppose that the assumptions of Theorem \3.S\ are satisfied and fix T > 0. As- 
sume moreover that T^ t = {7^ t}fe>i € (Bk>i ^(lP> ^ s a point of the sequence T^j = 
i^NVk^l w tth respect to the topology T pro< i. Then 

Tr(l-A 1 )...(l-A k )^] t <C k (7.28) 

for allk>\ and t G [0,T]. 

Proof. Theorem 17.31 follows from the higher order energy estimates of Theorem 17. 1( the proof of this 
fact can be found in |12| Proposition 6.3]. □ 



8 Convergence to the infinite hierarchy. 

In order to prove Theorem 13.21 we need to prove the convergence of the BBGKY hierarchy towards 
a hierarchy of infinitely many equations. In the argument, we will make use of the apriori bounds 
from Propositions 15.31 and Theorem 17.31 for k = 2. 

Theorem 8.1. Suppose that the assumptions of Theorem \3.2\ are satisfied and fix T > 0. Suppose 
that = {7^^}fc>i £ 0fc>i C{[0,T],£l) is a limit point ofT N}t = {j^ t } k=1 with respect to the 
topology T pr0 (i. Then T^t is a solution to the infinite hierarchy 



7S = W W (i) 7 LV 8 ^E/ dsU^(t-s)Tr k+l 5( Xj -x k+1 ) n ^ 



with initial data 7^0 = \(p}(tp\® k . Here lA^ k \t) denotes the free evolution of k particles defined 
in j^.5| ). 

Proof. Fix k > 1. Passing to an appropriate subsequence, we can assume that, for every 



sup TrjW 

te[o,T] v 







as N — > 00 . 



(8.2) 



We will prove (|8.ip by testing the limit point against a certain class of observables, that is dense 
in fC k . To characterize the class of observables we are going to consider, we define, for an arbitrary 
integer k > 1, 



We will consider £ IC k such that 



HS 



< OO, 



(8.3) 



where ||^4||hs denotes the Hilbert-Schmidt norm of the operator A. Note that the set of observ- 
ables j( k ' satisfying the condition (|8.3p is a dense subset of IC k . Moreover, using the fact that 



e l ^ t ( X j)e' 



-iA,-t 



( X j — itVj), it follows that 



n 7 k u {k \t)j^n 7 k <c{i + \t\) 7 n 7 k j^n 7 k 

HS 



HS 



(8.4) 
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Note also that, with the norm I J W I defined in (|6.3p . we have 

<C k 



HS 



(8.5) 



for a constant C k only depending on k (see [12] . Eq. (7.8)). Combining (|8.4j) with (|8.5j) . we also have 



\\u^{t)j^\i <c k {i + \t\) 7 n 7 k j^n 



HS 



(8.6) 



In order to prove Theorem 18.11 it is enough to show that, for every JW g £ fe satisfying (|8.3p . 

TrJ (fc) 7So = TrJ (fc) |^)(^r fc (8.7) 



and 



Tr J« 7 g t = Tr J^( fc )(t) 7 g - 8 W £ / dsTr J^ k \t - s) \s( Xj - x k+1 ),^ 



for all t G [0,T]. 

The relation (|8.7p follows from the assumption (|3.5p and from (|8.2 



In order to prove (|8,8p . we fix t £ [0, T], we rewrite the BBGKY hierarchy (|2,ip in integral form 
and we test it against the observable j( k \ We obtain 



Tr JW 7 « =Tr j( fe )w«(i) 7 ^ -iV / ds Tr J« ^ fc )(t - s)[Vj 
-i(N-k)J2f dsTiJ^U^ k \t - s)[Vi 

„-_i jo 



■9) 



jv^j -a?fc+l),72j ] 



From (|8.2p it follows immediately that 



Tr jW 7 « ^Tr J^j, 



(fc) 
oo.i 



5.10) 



and also that, as iV — > oo, 
Tr jW^W^hg^Tr 



-i 



(fc) 



ill 



- Tr (w«(-t)J«) 7 W = Tr J« W«(th£i, . 

Here we used that, if j( fc ) G /C fc , then also U {k \-t)J^ G £ A . 

Next we consider the second term on the r.h.s. of (|8.9p and we prove that it converges to zero, 
as N — > oo. To this end, we note that, setting = U^ k \t)J^ k \ we have 



£ Tr jW^( fc )(t - ,)[^(xi - Xj),^] = £ Tr jW [^(x* - ^), 7 Jgj 

and therefore, from (|6.12p and (|8.6p . we obtain that 

fc 



(fc). 



^ Tr jW^W^-s)^^-^),^ 



< Cfc III jjS III < C fc (l + T 7 )||^jW^|| H s 



i<j 



N 



N 
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for all < s < t < T. This implies that for all A; G N, for all t G [0, T] and for all jW such that (1831) 
is true, 



i<j 



Y, I dsTrJ^U^(t-s) VNixi-Xj),^ 
7?i Jo 1 







1.12) 



as N — > oo. 



Finally we consider the last term on the r.h.s. of (|8.9|) . First of all, we observe that, for every 
fc G N, t G [0, T] and j( fc ) G /C* such that is satisfied, we have 



kV] / d S Trj(%«(t- S )[v A r(x i -x fe+1 ),7^ 1) 

3=1 J ° 1 

as Af — > oo. This follows (similarly to (|8.12p ) from (|6.17p . and from the bound (|8.6p . 

It remains to show that, for every fixed k G N, t G [0, T], and for every G /Cfc with (|8.3 

iVV / daTrjWu^^t-sftVNfa-Xk+i),^] 
Tji Jo 



1.13) 



3=1 



k t 

SvraoV / dsTr J^U^(t - s)[5( Xj - x k+1 ),^] 

3=1 J ° 



5.14) 



as iV — > oo. To prove (|8.14p . we fix s G [0, t], and we consider, for example, the contribution with 
j = 1. We write 



Tr J( fc ^( fc )(t - sJJVVHsa - Sfc+OT^ 



Tr jf} t NV N {x x - x k+ i)W N ^ k+1) W* N ^ k+1) ^ 



(fc+i) 



(8.15) 



where W'jv'Ji.ik+l) denotes the wave operator associated with the Hamiltonian —A + (1/2) Vjy acting 
only on the variable x k+ \ — x\ (as defined in (|5.2p ). Therefore, if we choose a probability density 
h G L 1 (IR 3 ), with h > 0, J dxh(x) = 1, and we denote h a (x) = a~^h{x / a) for all a > 0, we have 

Tr jW\t - s)NV N ( Xl - x k+1 )^ + s 1] - SvraoTr J^ k \t - a)S(x 1 - x k+1 )^+^ 



< 



Tr J { ^\ [NV N {xi - Xk +i)W N ^ k+l) - 8vra 5{ Xl - x k+1 )] W* N ^ k+1) -f 



(fc+i) 

N.s 



+ 87rao 
+ 87rao 
+ 87rao 
+ 87rao 



Tr 4 k J t [S(x! - x k+1 ) - h a {x x - x k+ i)] W^ )(1)fc+1) 7^' 



(fc+i) 



Tr 



4-t h <*( x i ~ x k+i) W,(i,fe + i) - !)t]v," 



(fc+i) 



1.16) 



Tr jf\ h a (xi - x k+l ) (7^ 
Tr [h a (xi - x k+1 ) -S(x 1 - x k+ x)} 7^+ x) 



ffc+l) 



7, 



(fc+i) 



Here we insert the wave operator Wjvnfc-KH, because we only have a-priori bounds on the quantity 

W N,(l,k+l)' y N^ • Then we replace NV N (xi - cc*:+i)WjV,(i,fc+l) by 8Tra 5( Xl - x k +i). Afterwards, 

in order to remove the inverse wave operator ^ and to take the limit — ► 75^ , we 
need to replace the delta-function by the bounded potential h a independent of N. At the end, h a is 
changed back to the 5-function. 
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In Lemma 18.21 and Lemma 18.31 we prove that, for every k G N, for every < s < t < T, and for 
every j( fc ) G K h with (1831) . 



Tr J s ( 5 [JVVKr(a?i - »fe+l)^Ar,(i,fc+l) - 87ra 6(2:1 - x k +i)] W^,(i, fc+ i)7 
as A" —* oo and that 



(fc+l) 



as a — > 0, uniformly in iV. 

As for the third term on the r.h.s. of (|8.16p we remark that, for fixed k & N, s E [0, T], 6 /C^ 



(fc+l) 



.17) 
.18) 



and a > 0, 



Tr JS ^(xx - x fc+ i)(W^ i(lifc+1) - l) 7j (fc+ ' ' 



IN,s 



1.19) 



as N —>■ oo. In fact, for the bounded operator A = Jg_ t h a (xi — x k+ \), we can use the spectral 



decomposition 7^ = Ylj ^j\€j K '~ rXJ ){€j K ~ rX} I wrtn Aj = 1, Xj > 0, = 1, and estimate 



;(fc+l)wt(fc+l) 



(fc+l) 



^^(^ i(life+1) -l)7^ 1) < II^IIE^IIW.CW)- 1 )^ 



(fc+l)i|2 



„( fc + 1 ) 



<C\\A\\N- 1 / 3 Tr(l-A 1 -A k+1 ) 7kii 
< C\\A\\N 2 / 3 {i{j n , s , (H N + N)ip NiS ) < C||A||JV-V3 



(8.20) 



by the energy conservation and (|3,ip . From the first to the second line we used Lemma I8.51 Since 
the operator A is bounded for any fixed and a > 0, we obtain (|8.19p . 

To control the fourth term on the r.h.s. of (|8.16p we observe that, for arbitrary 5 > 0, 



Tr4%h a ( Xl -x k+1 ) ( 7 S: 1) -7^ 1) 
= Tr J;%h a (xi ~ x k+ i) 



1 + 6(1 - A fc+ i)V2 



(fc+i) _ (fc+i) 

IN,s "°o,s 



1.21) 



+ Tr jf\ h a {xi - x k+x ) ( 1 



l + «5(l-A fc+1 )V2 



Ak+1) _ Ak+l) 
iN.s loo.s 



The first term on the r.h.s. of the last equation converges to zero, as A" —* oo, for every fixed 6, a > 0; 
this follows from the assumption (|8.2p and from the observation that J^\h a (x\ — x k+ i)(l + 5(1 — 
Afc+i))" 1 is a compact operator on L 2 (R 3 ( fe+1 )). As for the second term, we notice that it can be 
bounded by 



Tr jf\ h a (xi - x k+ x) 1 



(fc+l) _ ( fc+ i) 

IN.s loo.s 



l + 5(l-A fc+1 )V2) 

< «||J«|| ||^||ooTr|(l - A W )^ ( 7 ^) - 7 (5+D) 

< 5a- 3 \\jW\\\\h\\ oc (Tr (1 - A fe+1 ) 7 ^ 1} + Tr (1 - A k+1 ) 7 ^ 

< C5a- 3 



1.22) 



uniformly in N. Choosing, for example, 6 = a 4 , it follows that 



Tr h a (xi ~ x k+ i) ( 7)^ 



_ (fc+i) 

IOO.S 



< r/(a,AT) + Ca 



1.23) 
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where r/(a,N) — > as N — > do, for every fixed a > 0, and where the constant C only depends on 
J( fc ). 



Finally, using Lemma [10,3l and Theorem 17.3^ the last term on the r.h.s. of (|8.16|) can be controlled 



by 



TrjS [/>„(,•, l )-,$(.,-, -, -,,+!)] ^ + < Cce 1 / 2 I Jni||Tr(l-Ai)(l-A fe+1 

<C(k,T, j( fc ))aV 2 . 



(fc+i) 



»V2 in 



i) 



5.24) 



From dSHSJ), (HUD, ([EgD , (HUD, ([833]) . and (|Q3D it follows that, for every jfe > 1, < s < i < 
T, and j( fc ) € /C fc with (IQ1) . 



Tr J^ k \t - s)NV N ( Xl - x fc+1 ) 7 S+ 1) - SvmoTr J«Z/ fe )(i - a )«5(xi - x k+1 ) 7 £f ] \ - (8.25) 
as N — > oo. Similarly to (|8.25|) . we can also prove that 

Tr J( fc )^( fc )(t - S ) 7 K 1) W JV (x 1 - x fe+1 ) - 8™ Tr J^ k \t - sftg+^Sfa - x k+l ) -» (8.26) 



as TV — > do. Since (|8.25p and (|8.26p remain valid if we replace x± by any Xj, j = 2, . . . , k in the 
potentials, it follows that, for every k > 1, < s < t < T, and G /C fe with l|Ol . 



^ (Tr j(%«(i - s ) [jV^i - x k+1 ),^ 

- SvraoTr j( fc )^ fc )(t - s) \s( Xj - x k+1 ),^ 



1.27) 







asA^^oo. FromdEIID (with j( fc ) replaced by U {k \s-t)J^ k \ using the fact that ||W( fc )(5-t)j( fe )| < 
Ct for all < s < t < T), and from an estimate similar to (|6.17p but with "f^^ replaced by 760^ 
and NVpjixj — x k+ \) replaced by 6(xj — x k+ %), we can now apply the dominated convergence theorem 
to conclude <KW\ . □ 



The following lemmas are important ingredients in the proof of Theorem 18. II 

Lemma 8.2. Under the same assumptions of Theorem \8.1l and using the notation jf } = U (k \t)J (k \ 
we have, for every k > 1, i = 1, . . . , jfe, < s < t < T, J( fc ) G K k such that (ED is satisfied, 



Tr j(% [NV N {x e - x k+ i)W N ^ k+l) - 8vra S(x e - x k+1 )) W* N ^ k+l) ^ 



(fc+i) 

s 



as N — ► 00. 



Proof. We fix •£ = 1. Decomposing 7^^ 1 ' 1 = -\? 
u = (x± + Xfc + i)/2 and v = x\ — Xfc+i, we find 



t (fc+i) wt (fe+i). 



and introducing the variables 



Tr J s (fc) t NV N (xi - x fc+ i)7^ 1) = ^2>^j J dudvdx 2 ■ . . dx fc dx' fc J s ( ^(x' fc ; u + v/2, x 2 , . . . , x k )NV N (v) 



x 6 k+1 \u + «/2,x 2 , . . . ,x k ,u- v/2)?r i) (y! kt u- v/2). 



(fc+i), 



3.28) 
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The potential Vn(v) forces v to be of order 1/N. Using this fact, we are going to remove the 

— (fc+i) 

^-dependence from the observable and from the wave function £j . After removing this in- 
dependence, we introduce the wave operator using its L 2 -unitarity. We find 

Tr 4 k J t NV N (x 1 -x k+1 )^+ 1) 

= Aj J dudvdx 2 ..dx k dx' k NVn(v) ^ k+1 \u + v/2, x 2 , ..,x k ,u — v/2) 

3 

X [ J s ( -k x fe; u + v l% x 2> ■ ■ ■ j x k)lf +1) (x' fc , u - v/2) - jj^x'j.; u,x 2 ,..., x k )^ +1 (x* 
+ y~l Aj / dududx 2 ..dx fc dx / fc iV(W^yAr)(t;) J s ( ^(x' fe ; it, x 2 , ..,x k )£ k +1) (x' fc ,u) 

X [( W/ Ar,(l,fc+l)^' C+1) )( U + W A X 2, ••, X k ,U- V/2) - {W* N ^ k+l) if +1) ){u, X 2 , ..,X k ,u) 

+ Yl ^ {^j dv ( W N V n)(v)J J dudx 2 ■ ■ ■ dzfcdx'^ jj® t (x! k ; u,x 2 ,...,x k ) 
x if +1) (4> u )(^,(i,fc+i)^ +1) )(u,x 2 , ...,x k ,u). 

(8.29) 



From Lemma 18.41 we know that 

J dvN(W^V N )(v) = 8™ . 

Therefore, from (|8.29p . we obtain that 

Tr jf\ [NV N (xi - x k +i)W N> ( 1}k+1 ) - 87ra (5(xi - x k+1 )] W^ ljk+1 ^^ 1 ' ) 

< J dudvdx 2 . . .dx k dx' k NV N (v)\^ k+1 \u + v/2,x 2 , . . . ,x k ,u - v/2)\ 

j 

x pi-tC**; u + v / 2 > x 2 ,..., x k )lf +V> (x! k , u - v/2) - J s (fe) t (x' fc ; u, x 2 ,..., x k )^f +1 \x! k ,u) 
+ y^Aj / dud-udx 2 • • • dx fc dx' fc N(W^V N )(v) (x' fc ; u,x 2 ,..., x k )^f +1 \x k , u) 

3 J 

X [W,(l,fc+l)Cj fc+1) )(« + U A S2, • • • , Xfc, « - V/2) - (W^ )(lifc+1) ff +1) )(U, X 2 , • • • , Xfc, u) 

=: ^A ; (I, • ) . 

3 

(8.30) 

The terms L can be bounded exactly like the term I on the r.h.s. of (16.141) . after replacing by 
J {k \. Following the steps (lUT^D - dUT^D . we find that 

C IIU (fc) lll 

Y, X i l 3 < Tr (" A i " A ^+i + !) ^(i^+DT^^^d.fe+i) (8-31) 

j 

which converges to zero as iV — > oo for all < s < t < T, and all observables satisfying (I8.3() 
(here we used the a- priori estimate given in Proposition 15.31 and the observation (|8.6p ). 
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Next, we consider the second term on the r.h.s. of ()8.3Qj) : 



II; 



dudvdx 2 ■ . . dx k dx' k jj® t (x' fe ; u, x 2 , ■ ■ ■ , Xk)^^ ( x fc> 

x (N 3 (W*V)(Nv) - S(v)) (W^ k+1) ^ k+1) )(u + v/2, x 2 ,...,x k ,u- v/2) 



5.32) 



To control this contribution, we first insert a cutoff x( v )i this wm allow us to apply Lemma 110.21 
to bound the integral over u and v. To this end, we choose a function x S C^°(IR 3 ) such that 
< xi x ) — 1; x( x ) = 1 f° r \ x \ — 1 and x( x ) = for M — 2, and we put X = ^ ~ X- Using Xi we 
decompose the r.h.s. of (|8.32|) in two parts 



II; < 



dudvdx 2 . . . dx k dx.' k J s ( ^(x' fc ; u,x 2 ,... , x k )x{v)£^'' {* k , u) 

x [N 3 (W*V)(Nv)-6(v)] (W* N ^ k+1) ^ k+1) )(u + v/2,x 2 ,...,x k ,u-v/2) 



=(*+!) , 



+ 



dudvdx 2 . . . dx fc dx' fc J {k \(^ k ;u, x 2 ,..., x fc )x(w)^ _r±; (x' i ., u) 

x N 3 (W*V)(Nv)(W* N ^ k+1) ^ k+1) )(u + v/2, x 2 ,...,x k ,u- v/2) 

=:Aj + Bj. 

The term Bj can be bounded by 
Bj ,< J dudvdx 2 . . . dx fc dx' fc jS(x' fc ; u, x 2 , ■ ■ ■ , x k ) X (v)N 3 (W*V)(Nv) 

x fkf +1) (x' fc ,n)| 2 + (W* N ^ k+1) ^ k+1) )(u + v/2,x 2 ,...,x k ,u-v/2) 



=(*+!), 



(8.33) 



< lie 



(fc+l) N 2 



sup / dx 2 • • • dxfc 


r(*0 / / 











|d|>JV 



|(W*V)(u)|dt; 



+ (sup y dx' fc ^(x'^Xfe) 

x J dxi . . . dx fc+1 x(xi - xjfc+OiV^WV)^^! - x fc+ i))| (^ >(lifc+1) ^ (fc+1) )(x fc , Xfc+i) 
From Lemma llO.lt we obtain 



j \J\v\>N 



dv \{W*V){v)\ 



< (Tr ((Vi • V 2 ) 2 — Ai — A 2 + l) W N ^ k+1) 4+ l) W* N ^ k+l) 



5.34) 



as iV ^ oo. Here we used Proposition 15.31 and the fact that, since W*V E L 1 

[ \W*V(x)\dx -» as iV^cx). 

J\x\>N 



As for the term Aj on the r.h.s. of (|8.33p . Lemma 110.21 implies that there exists a sequence 5n — * 
as A" — > oo (5n corresponds to the sequence 0\m defined in Lemma [W.2\ with V replaced by W*V) 
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such that 



dudv 



,1/2 



Aj < Sn J dx 2 . . . dxfcdx' fc 

((A u - A„) 2 - A u - A v + 1 

(1 - A u + A 2 v ) 1/2 x(v) Jst\(x' k ; u,x 2 ,..., x k )lf +1) (x' fc , u) 
j dx 2 ...dx k dx k (| du|(l-A u ) 1 /2|( fc+1 )( x ' fc)U )| 2 y 2 
J s (A : ) 4 (x / A .;'u,X2, • • • ,x fe ) + V u J^ t {x k ;u,x 2 ,...,x k ) 



2 \ 1/2 



dudv 



< SnWxWep 
x sup 



x ( / dudv 



{(A u - A v f - A u - A v + 1) 1/2 

x (^7V,(i,fc+i)£f + V A sea, • • • . x k , u - v/2) ) 



2 \ I/ 2 

(8.35) 



With a Schwarz inequality, we find 
Aj < 8 n \\x\\h* J dx 2 . . . dx fc dx' fc 



x sup 



J^} t (x' k ;u,x 2 , • • • ,x k ) + V u j(%(x' k ;u,x 2 ,...,x k ) 
Qdu\(l-A u )^f +1 \ X ' k ,u) 2 



+ J dudv\ {(A u - A v ) 2 -A u -A v + 1) 1/2 (W^ )(1)fe+1) ef (u + v/2,x 2 , ..,x k ,u- v/2)Q 



<C(k,T,J^)6 N ((^* +1, ,(l-A B )ej* +1) ) 
^v,(i,fc+i)£j 



+(n.n. fc+ i^f (cvi • v fc+1 ) 2 - Ai - A fe+1 + 1) n, ( i, fc+ D# +1) ; 



From Lemma 110.21 and Proposition 15.31 we find 

^ A, Aj < C^at as N -> oo 



(8.36) 
(8.37) 



and this, with (|8.34|) . implies that 



V] AjIIj as iV -> oo . 



Together with (|8.3ip and (|8.30p , this concludes the proof of the lemma. 



□ 



Lemma 8.3. Under the same conditions as in Theorem \8.1\. we have, for every k > 1, i = 1, . . . , k, 
< s < T, and J (fc) G K k satisfying (8.3\) . 



Tr 4-t [S(x t - x k+ i) - h a (xi - x fc+ i)] W^ i( ^ )fe+1) 7^' 



(fc+i) 
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as a 



0, uniformly in N. Here we use notation J t (fe) =U ( - k \t)J ( - k h 
Proof. We fix £ = 1. Using the decomposition = J2j Aj|£j fc+1 ^)(<!;j' C+1 ' ) |, we find that 

|Tr J s ( 5 [5{xi - x k+1 ) - h a (xt - xfc+i)] Wjv^i.fc+ijT^J" 15 

< £ | / ■ ■ ■ W4. «(4; » + .A *.,....*> HW - *.<•>] 

x ?f +1) (x' fe , u - v/2)(W* N ^ k+1) if +l) ){u + v/2, x 2 ,...,x k ,u- v/2) 

< ^2 | / dudvdx2 ■ ■ ■ dxfcdx' fc h a (v) 

X [ J S( x fc; u + W A x 2 ,---, x k )tf +1 \x' k ,u - v/2) 

x + V A ^2, ■ ■ ■ , x k ,u - v/2) 

- J { s k J t (x' k ;u,x 2 , . . . ,x k )tf +1 \x k ,u)(W^^ hk+1) ^ k+1) )(u,x 2 , ...,x k ,u) 



Similarly to (|8.29[) . we first replace v by in J^} t (x! k ;u + v /2,X2, . . . , x k )£,y~ r '~' (x.' k ,u — v/2) and then 



(fc+i) , , 



in (W* NAhk+1) ^ k+1) )(u + v/2, x 2 ,...,x k ,u- v/2). 



We obtain 



Itt 4 k _\ [<5(xi - xfc+o - /» tt (xi - x fc+ i)] - " 



JV,(l,fc+l)7iV,s 

< ^Aj| / dndudx 2 ...dx fc dx / fc /i Q (i')(^ (liA . +1) ^' c+1) )(n + u/2,X2,...,a;fc,u-v/2) 
x [jS(4; u + x 2 , . . . , x k )^f +1 \x' k , u - v/2) 



-ji%(x' k ;u, x 2 ,..., x k )^f +1 \x' k ,u) 



+ y~^Aj I / dnd-udx 2 . . . dx fc dx' fc 7i a (u) J^ t (x' k ; u,x 2 ,..., x k )cf +1 \x' k ,u) 
j J 

x [( W N,(l,k+T)£j k+t) )( u + u / 2 ' ^2, • • • , x k ,u - v/2) - (Wx : ( l!k+1) £j h+1) )(u, x 2 ,..., x k ,u) 



(8.38) 



To bound the first term, we expand the difference in an integral 

Js-tftki u + v / 2 > x 2, Xk)if +l) {* k , u - v/2) - J g ^(x4;u, x 2 , ■ ■ ■ , Xk)^^ 1 (x' fc , «) 

.1/2 



^ dr w • Vi J s (fc) t (x' fc ; u + rv, x 2 , . . . , x k )^f +1) (x' k , u - rv) 



(8.39) 



1/2 



dr J { s k J t (x' k ;u + rv,x 2 ,..., x k )v ■ V k+1 ^ ±J (x' k ,u - rv) 



(k+l), , 
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and we obtain that 

IIIj < J dudvdx 2 ■ ■ ■ dxfcdx'fc J dr h a (v)\v\ {W^j k+1 ^)(u + v/2, x 2 , ■ ■ ■ , Xk, u - v/2) 
x ( ViJ s (fc) t (x / A .;'u + r?;,X2,...,x fe ) \^ k+1) (x' k ,u - rv 

+ Js-t( x 'k'i u + rv,x 2 ,..., x k ) V k+1 ^ k+1 \x' k ,u - rv) ) 

which implies that 

III, < a C(k, T, J«) (Tr (1 - A fc+a ) 7 ^ 1) + Tr ((Vi • V fc+1 ) 2 - Ax - A fc+1 + l) 7 <*+ 1 > 



.40) 



<C a. 



The terms IV,- can be estimated similarly to the terms II, considered in ()8.32|) ; in particular, analo- 
gously to (grjEg) and ([5757]) . we also find 



2 A i IV i < % T > jW ) Pa (iV (1 - A fe+1 ) 7 {f+ 1) + Tr ((Vi • V fc+1 ) 2 - At - A k+1 + l) 7 # 



(fc+i) 



<C/3 Q 



where /3 a — ► as a — > uniformly in N (the sequence (3 a comes from Lemma 110.21 with V replaced 
by h). This concludes the proof of the lemma. □ 



Lemma 8.4. Suppose that V > 0, with V(x) < C(x) a for some a > 5 (this implies, in particular 
that V <E L 1 ^ 3 ) n L 2 (R 3 ) and thus that V is in the Rollnik class of potentials). Let W denote the 
wave operator (as defined in Proposition \5.1\) associated with the Hamiltonian f) = — A + (l/2)V(x). 
Then 

r dx (W*V){x) = 8to 



where oq is the scattering length of the potential V . 

Proof. First of all, we observe that, under the assumption that V > and V(x) < C{x)~ a , for some 
a > 5, the operator f) = —A + (1/2)1/ cannot have a zero energy resonance (recall that a zero-energy 
resonance of f) is a solution </> of (—A + (l/2)V)ip = such that |v?(^)| < C/\x\ for all x G M 3 ); this 
can be proven using the maximum principle. We will make use of this observation in the proof of 
this lemma. 

Next, we note that, since W* maps L 1 (M 3 ) into L X (]R 3 ) (see Proposition 15. lh . we have that 
(W*V) e L 1 ^ 3 ) and thus 



dx (W*V)(x) = lim / dx (W*V)(x) X e(x) = lim / dx V(x)(W X s)(x) 



.41) 



with 



Xe{x) 



1 



1 + ex 2 

We expand W\e m terms of solutions <p(x, k) of the Lippman-Schwinger equation 

1 



<p{x,k) 



Ak-x 



e i\k\\x-y\ 

. , d y-, rV(y)ip(y,k). 

8ir .1 \x - y\ 



.42) 
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It follows from |24| Theorem XI. 41, a)] that Eq. (|8.42p has a unique solution (p(x,k), such that 
f(x, k)V l l 2 (x) G L 2 (M 3 ), for all k G M. 3 such that k 2 g" £, for an exceptional set £ with Lebesgue 
measure zero. The set £ consists of all values of k 2 for which zero is an eigenvalue of the operator 



3.43) 



From the observation that the operator f) = —A + (1/2)1/ does not have a zero energy resonance, it 
follows immediately that G" £; in fact, if Mqi/j = for some ip G L 2 (M 3 ), then 



■0(x) 



V/2 (5 



f - y| 



which implies that VK^)/^ 1 ^ 2 ^) < f° r \ x \ ^ 1 an d thus that tp(x) := ip(x)/V 1 ^ 2 (x) is a zero- 

energy resonance solution of (—A + (1/2)V) ip = 0. Since Mo is a non-negative Fredholm operator 
with no eigenvalue at zero, it follows that there exists A > with <t{Mq) C (A, oo) (here cj{Mq) 
indicates the spectrum of Mo). Since moreover — Mq is a compact operator with kernel 



1 % I I *C y | 1 

{M lkl -M ) (x;y) = -V^ix f ~ V l ' 2 (y) 



f - y| 



we obtain that 



I Mi 



1*1 



M 



OllHS 



dxdy V(x)V(y) 



,i\k\\x-y\ _ j| 

|x — y| 2 



< 



|£f||F» 2 



L 1 



3.44) 



3.45) 



and thus that there exists k > such that a(M\ k \) C (A/2, oo) for all \k\ < k. In particular it follows 
that 



IM, 



|fc| 



< 2/A 



for all fc£l with Ifcl < k. 



From [24l Theorem XI. 41, e)l we also find 



(Wxe)(x) = L.I.M. (2vr)- 3 / 2 / dk <p(x, k) %{k) = L.I.M. / dA; tp(x, 



k). 



-\k\/V^ 
&7r\k\e 



(8.46) 



3.47) 



< 



where L.I.M. denotes the L 2 -limit as M — > oo and 5 —> of the integral over {k G IR 3 : 
M and dist (k 2 ,£) > 5}. Inserting (|8.47p in the r.h.s. of (|8.4ip . we find (recalling that k > is 
chosen such that (18. 46ft holds true) 



/ 



dx (W*V)(x) Xe{x) 



e -\k\/V^ 

dxdk V(x)ip(x, k) — ; h 

\>k ' 47r|fe|e 



e -|*l/v^ 

dxdk V(x)ip(x, k) — - — 

|fc|<« 47r|fe|e 



3.48) 



The first term on the r.h.s. of (|8.48[) can be controlled by 

e -\k\/Ve 



/ dxdkV(x)(p(x, k)- 
J\k\> K ' 47r|/c|e 



dk V^k)'- 



>K 



4n\k\£ 



< C 



< c\\v\ 



I? 

'|fc|>K 

e -K/(2y/e) 



dk- 



-2\k\/s/i 



\ 2 E 2 



1/2 



3.49) 



L 2 ' 



£ 



3/4 



37 



as e — > 0. Here we introduced the function 

V\k) = l.i.m.(2^)" 3 / 2 / dxV(x)f(x,k) 



where l.i.m. denotes the L 2 -limit of the integral over |x| < M as M — > oo; the existence of V" for 
V G L 2 (M 3 ) and the fact that ||V'||^2 < ||V 11^,2 (actually, in our case, ||^"|| i2 = ||V[|j^2) are proven 
in [24\ Theorem IX. 41]. As for the second term on the r.h.s. of ()8.48p . we have 



e -l*IA/s 

dkdx V(x)(p(x, k) — — 

|*|<« A7r \k\e 



\k\<K 



dkdx V(x)(p(x, 0)- 



-\k\/V~e 

A-K\k\e 



e -\k\/y/i 

+ / dkdx V (x) (<p{x,k) - ip(x,0)) ——— 
\k\< K 4vr|fe|e 

1- {1 + K e- 1/2 )e- K£ ~ 1/2 ^j J dxV(x)<p(x,0) 

f e -\ k \/\/£ 

+ / dkdx V(x) {ip(x, k) - ip{x, 0)) 

J\k\<K 4Tr\k\e 

Using that (p(x, 0) is the solution of the zero energy scattering equation 

(-A + (l/2)U(x))v?(x,0) =0 

with the boundary condition (p(x,0) -> 1 as \x\ — > oo, it follows that (see (|1.6|) ) 

dx V(x)ip(x, 0) = 87rao • 



To bound the second term on the r.h.s. of (|8.50p . we define 

Mx) =V 1/2 (x)ip(x,k) 
and we observe that, from the Lippman-Schwinger equation (|8.42p . 



^ k {x)-Mx) = V l, \x){e ik - x - I 



irhA^LAvWiy) (Mv) -My)) 



[ ayv 1/2 (x) e - 

8vr J \x-y 



-I / dyV^( x y- 
8vr J \x-y\ 



v 1/2 (y)My), 



which implies that (with Miju defined in (|8.43|> ) 



[M| fc , - M] (x) = V 1/2 {x) (e lk - x - l) - i- / d//T ' 



i|fc||x — y\ i 

/ 2 (z) , V 1 /- 



f - v\ 



By dSHBJ), we have 



L 2 



< c 



V^fa:) (e** - 1 



L 2 



+ \\v 1/2 \\ L i\k\ / ^ 1/2 (y)l^o(y)| 



<C\k\(\\\x\*V\\]!U\\V\\f) 
< C\k\ . 



1.50) 



5.51) 



.52) 



(y)My)- (8-53) 



.54) 
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Therefore, the second term on the r.h.s. of (|8.50p can be bounded by 

e -[k[A/E 



dkdxV(x) (<p(x,k) - <p(x,0)) 
\k\< K 4ir\k\e 



< / dkdxV x / 2 (x) \ipk(x) — fpo{x)\ . ... 



<||V||K 2 / dfc||^-Vo||L 2 : 



-|fc|A/i 



< CM 1 * 

< Ce 1 ' 2 



\k\<K 47r|A;|e 
e -|*l/VE 



dA:- 

|fc|<K 



5.55) 



and thus it converges to zero as e — ► 0. The last equation, together with (|8.4ip . (|8.48[> . (|8.49[> . (|8.50[> . 
and (18.5ip . concludes the proof of the lemma. □ 

Lemma 8.5. Suppose that V > and V{x) < C{x)~ a , for some a > 5. Then, for every g £ 
L 2 (R 3 ,dx), we have 

"{W N -l)g\\ < CN-^Mh! . 
Proof. Let t) N = -A + (l/2)V N (x). Since 

W N = s - lim e ii)Nt e lM 



t— >oo 



it is enough to prove that 



sup 



e -it, N t _ e iAt g 



Note that 



which implies that 



d 
dt 



e -it, N t _ e iAt g 



< CN- 1/6 \ 



21m (e-^g^Nixy^g) 



9\\m 



1.56) 



iAt. 



< 2 / ds 



(e- i ^ s g,V N (x)e lAs g) 



Next we observe that 

<e-^ s 5 ,%(x)e iA * 5 >| < He-^Hoo ||e iAs 5 ||oo Mi < 
where we used the fact that 



Ns 3 



3.57) 



3.58) 



IIWjvII 



p^p — 11^^ lip— *p 



for every N and 1 < p < oo. For small s we need a different estimate of the integrand on the r.h.s. 
of (|8.57p . To this end we remark that 

(e-^ s g,V N (x)e lAs g}\ < (e^ s g, V N {x)e l ^ s g) 1 ' 2 {e lAs g, V N (x)e lAs g) 

<C||^v||3/ 2 ||Ve^ s <7||||Ve iA ^|| (8-59) 

<c||y|| 3/2 (i + ||y||3/2) 1/2 lbll^, 
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where we used that 1 1 Vjv 1 1 3 /2 = IIVII3/2 and we estimated 

HVe^sH 2 = (e l ^ s g, -Ae^g) < (g, t> N g) < (1 + \\V\\ 3/2 )\\gf Hl . 
Combining (|838j) and (jH39|) . we obtain from (f8~57|) 



e ^Nt_ e iAt ]g 



< 



2/ d S ||F|| 3/2 (l + ||F||3/2) 1/2 

J 



\g\\m 



+ 2 1 ds 
In 



\V\\i\\g\\ 2 Hl ||W||o™ 



Ns 3 



< ( CiN~ a + C 2 N 2a ~ 1 



\g\?m 



for every t £ TSL. Choosing a = 1/3, we obtain (|8.56p . 



$.60) 



□ 



9 Approximation of the initial data 

In this section we show how to regularize the initial wave function V'jv given in Theorem 13.21 

Proposition 9.1. Suppose that ipN G L 2 (M. 3N ) with HVwII = 1 is a family of N -particle wave 
functions with 

(^ N ,H N ^ N ) <CN (9.1) 
and with one-particle marginal density 7!^ such that 

7$ - k><¥>| asiV^oo (9.2) 
for a f £ H 1 (M 3 ). For k > we define 

X {kH n /N)^ n 



\x(kH n /N)iI>n\ 



(9.3) 



1 



is a cutoff function such that < x < 1, x( s ) = 1 f or < s < 1 and x{ s ) = for 



Here x £ C| 

s > 2. We denote by j)^' , for k = 1, . . . , N, the marginal densities associated with i/jn- 
i) For every integer k > 1 we have 



(il) N ,H k N ^ N ) < 



2 k N k 



(9.4) 



ii) We have 



sup HVw - <MI <Ck 1 ' 2 

N 



in) For k > small enough and for every fixed k > 1 we have 



lim Tr 

N^oo 



. 



(9.5) 
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Proof. For the proof of part i) and ii) see \12\ Proposition 8.1]. To prove in), we begin by noticing 
(see (jl.2p ) that it is enough to show that 



lim Tr 

jV-*oo 



0. 



Moreover, since the limiting density is an orthogonal projection, trace-norm convergence is equivalent 
to weak* convergence. In other words, it is enough to prove that, for every compact operator 
€ /Ci and for every e > there exists iVo = No(J^,e) such that 



Tr JW [7$-\<p){<p\ 
t by ol 

L^R 3 ^- 1 )), with = 1 such that 



< £ 



(9.6) 



for N > Nq. To show (|9.6j) . we start by observing that, from (|9.2p . there exists a sequence ^jf ^ £ 



b N ~ H> ® Of I 



as — > oo . 



(9.7) 



This was proven by Alessandro Michelangeli in [23] using the following argument. Choose an or- 

thonormal basis {/$}i>i of L 2 (M 3 ) with f\ = (p. Choose also an orthonormal basis {gj}j>i of 
£2( K 3(iV-i))_ Then Qne can write 

(N) 



IpN = E°i' fi®9j 



and 



This implies that 



hM(Viv|= E ^^l/iX/i'l®^-)^-'!- 



and therefore, using (|9.2|) . that 



EI-S3 



(AT) 1 2 



as iV — > oo. Thus, putting ^f ^ 



2Zj a i,j 9j, we get 



V - 9? ® £ 



2tJV-l)„2 



A' 



EEM?'i 2 = i-Ei< 



(JV),2 



as iV — > oo. It is then simple to check that £ 



(JV-1) _ *(JV— 1) /uTiN-l), 



N 



N 



satisfies ([97 



On the other hand, there exists tp* € i? (R ) with ||y>*|| = 1 and such that 

£ 



If — <P*\\ < 



32||JW|| ' 



Let E = x(kH n /N). Then 



H - 1)iPn\\ 2 < jr{ipN,H N ip N ) < Ck 
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independently of N. Therefore, choosing k > so small that ||H?/>jv|| > 1/2, we find 



A? 



< 

< 4 

< 4 

< 



, „ JN-l) 



(9.8) 



+ 4\\ip - cp*\ 



6||JW|| 

for all N sufficiently large. Next we define the Hamiltonian 

N N 

Hn ■= -^Aj + ^2 VN(xi-Xj). 

j=2 Ki<j 



(9.9) 



Note that Hn acts only on the last N — 1 variables. We set H := /N). Then we claim that, 

if e > is small enough, 



<p*®C N 



< 



3||JW| 



(9.10) 



for N sufficiently large. The proof of (|9.1U|) can be found in [12, Proposition 8.1]. To get (|9.6p we 
define 



4>N ■-- 



3 (<£* Cjf 1 

|s(^®^ _1) ) 



where we used the fact that S acts only on the last N — 1 variables and the fact that \\<p*\\ = 1. 
Define 

Note that ipN is not symmetric in all variables, but it is symmetric in the last N — 1 variables. 



In particular, 7^ is a density matrix and clearly j^' = \<p*){tp*\. Therefore, since \\tf) 
e/(3||jW||) by (I9TTU]) and since ||y>- ¥>*|| < e/(32|| J (1) ||), we have 

Tr J« (7$ - |^|) I < |Tr J« (7$ - |^)<¥>*l) I + |Tr J« (|^)(^| - |^|) 

< 2|| J {1) || H^iv - -0iv || + 2|| J (1) || H99 - <^*|| < e 
for iV sufficiently large (for arbitrary n, e > small enough). This proves (|9.6p . 



(i) 



A" 



(9.11) 



□ 



10 Poincare-Sobolev type inequalities 

In the proof of the convergence we need to estimate potentials converging to a delta functions, and 
their difference to a normalized delta-function. To this end we make use of the following three 
lemmas. 
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Lemma 10.1. Suppose V G L^IR 3 ). Then 

\(<p, V( X1 - x 2 )^)\ < C\\V\U (V, ((V! • V 2 ) 2 - A x - A 2 + 1) V) 1/2 Q 

x ( 9 9,((V 1 -V 2 ) 2 -A 1 -A 2 + 1) V 9) 1 / 2 

for every ip,ip G L 2 (M 6 , da?i, da^). 

Proof. Switching to Fourier space, we find 

(ip,V(xi - x 2 )ip) = J dp 1 dp 2 dq 1 dq 2 fi(pi,P2)'4>(qi,q2)V(qi-Pi)5(pi+p2-qi-q2)- (10.2) 
Therefore, by a weighted Schwarz inequality, 
(ip, V(xi - x 2 )ip) 

< Halloo ( [ dp 1 dp 2 dq 1 dq 2 ^ t ^2 t ^2 t ! » P2) | 2 5(pi +P2 - ?1 - 92)) 

V Wl • g2j 2 + 9i + <7 2 + 1 / 

X ( / dp ld p2d gl d g2 ^ • g2 f 2 I 1 tit! l£fo> ^)! 2 ^ +^ 2 - * - ^ 2 )) l/2 (10 - 3) 

vy (pi -P2) 2 +h +p 2 + 1 / 

sri11 (t/' V^-^ + ^-^-m ) 

x (ft ((V, • V 2 ) 2 - A t - A 2 + 1) </>) 1/2 (<p, ((V t ■ V 2 ) 2 - Ai - A 2 + 1) v ) 1/2 . 
The lemma will then follow from 

^ q {q ■ (p - q)) 2 + q 2 + (p - q) 2 + i < ^° ^ 



dq 



To prove (|10.4p . we proceed as follows. 

I M (q-{p-q)) 2 + q 2 + {p-q) 2 + l = L-i 



{q-ip- q)) 2 + q* + {p-qY + l J^ >w n _ g) 2 _ £\ 2 + g2 + {p _ g)2 + 1 



+ dg 



k-f|<bl (( g _2)^_2fj +(7 2 + (p _ g) 2 + 1 

(10.5) 



2 



The first term on the r.h.s. of the last equation is bounded by 

1 



/ d( ? rr~ ^ / d( ? 



l>H ((q-lf-f) +q 2 + (p-q) 2 + l 4-fl>bl ^ k — f| +1 



/ 16 /" j 1 

< — / do — -. < 00 



(10.6) 
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uniformly in p € R 3 . As for the second term on the r.h.s. of (|10.5p . we observe that 
/ dq 2 

i\<\p\ f^.E) 2 .^ +q 2 + {p _ q) 2 + l 



dx 



2 +(x + f) 2 + (x-f) 2 + l 



xr 



bl r 2 



47r / dr- 



2 -M!) 2 + 2r2 + M^ + i 



r 

/•H/2 
< C|p| 2 / dr 
J-M/2 



(10.7) 



bl/2 r 2( r + \ p \f + ( r + \fj +M! + i 

/■IpI/2 i /" 1 

<C dr^ -<C dr^ < oo, 

J-\p\/2 r 2 + l J R 



uniformly in p. □ 

Lemma 10.2. Suppose V £ L 1 (M 3 ) with J V(x)dx = 1. For a > , let V a (x) = a~ 3 V(x/a). Then 
there exists a sequence (3 a with (5 a — * as a — > such that 

\(<p, (V a ( Xl - x 2 ) - 5( Xl - x 2 ))i>)\ < C(3 a ty, ((Vi • V 2 ) 2 - Aj - A 2 + 1) V) 1/2 

x ((Vi - V 2 ) 4 + (Vx + V 2 ) 2 + 1) vp) 1 / 2 , 

for all ip,tp £ L 2 (R 6 ). 

Proof. Switching to Fourier space we find 

(tp, (v a (x 1 - x 2 ) - 5(xi - x 2 j^ip^ =1 + 11, 

where we defined 

I = / dpidp 2 dqxdq 2 dxV(x)^(pi,p 2 ) ( e iax <P^~^) _ {\ q 2 )5{p 1 + p 2 -q x - q 2 ), 

J|x-(pi-gi)|<a-i/2 V / 

II = / dpidp2dgidg 2 dx V r (x)^(pi,p 2 ) ( e iax ^ pi ~ gi) - l) $(q u q 2 )5(pi +p 2 - gi - £2) • 

■/|a>ft>i-ai')l>ar 1 /2 V / 



(10.9) 

To bound the first term we use that \e lK — 1| < \k\, for k £ R, and we observe that 



iii < a i/2 i^iii / dp ld p 2dgl d 92 y^ 2 'r + spg + 1 |^ 2 )| 



• «2) 2 + ?i + q 2 

v V (?i ■ q2) 2 + qf + q'i + l ,7, , 

Vbi -P2| 4 + (pi+p 2 ) 2 + i 

(10.10) 
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With a Schwarz inequality, we obtain that 

|I| <a 1/2 ||V||l ( I d Pl dp 2 d qi dq 2 \P^t^^^^ \ip(p 1 ,P2)\ 2 S(p 1 +p 2 -q 1 -q 2 )" 
\J [Qi ■ I2Y + q{ + q 2 + 1 

1 /2 

x ( f dp 1 dp 2 dq 1 dq 2 Q2 \ + ^ + Q2 ^ 1 = \${qi, q 2 )\ 2 5{pi + p 2 ~ gi ~ 92) ) 

< a 1/2 ||y||i L, ((Vi - V 2 ) 4 + (Vi + V 2 ) 2 + 1) <p) ((Vi • V 2 ) 2 - Ai - A 2 + 1) ^ 



\ 1/2 / \ 1/2 

dq \ ( f dp \ 



X 1 ' \q~p\ A +P 2 + l [»SsJ (p-(q- P )) 2 +P 2 + (q-p) 2 + l 



From 

f dq f dq 

peR»i \q-p\ i +P 2 + l J \qJTl 

and (fTO^j) it follows that 

|I| < Ca 1 / 2 ^, ((Vi - V 2 ) 4 + (Vi + V 2 ) 2 + 1) y?) 1/2 ((Vi • V 2 ) 2 - Ai - A 2 + l) ^ . 

(10.11) 

In order to control the second term in (|10.9p . we bound it by 

|II| < 2 / dpidp 2 dqidq 2 dx\V(x)\\(p(pi,p 2 )\\^(qi,q 2 )\5(pi+p 2 - qi- q 2 ) 

•V(pi-gi)l>"- 1/2 



<2 / dpidp 2 dq l dq 2 dx\V{x)\\(p{p 1 ,p 2 )\\'il){qi,q 2 )\5(pi+p 2 - qi- q 2 ) 

<\x\>a-y 4 



+ 2 / dpidp 2 dqidq 2 dx\V(x)\\(p(j) 1 ,p 2 )\\'4)(qi,q 2 )\5(pi+ p 2 - qi- q 2 ) 

J\pi-qi\>a--L/± 

< Pi,a J dp!dp 2 dqidq 2 \fi(pi,p 2 )\ \4>(qi,q2)\ $(pi + P2 ~ qi ~ 92) (10.12) 



+ 2||V||i / dpidp 2 dqidq 2 \fi(pi,p 2 )\\'4>(qi,q 2 )\5(p 1 +p 2 - q! - q 2 ) 

J\p 1 - qi \>a- 1 / i 

< p 1>a (if, ((Vi - V 2 ) 4 + (Vi + V 2 ) 2 + 1) ip) 1/2 (if,, ((Vi • V 2 ) 4 - Ai - A 2 + 1) V f 12 

+ 2||^||i / dpidp 2 dqidq 2 \(p{pi,p 2 )\\ : 4){q 1 ,q 2 )\5{pi+p 2 - qi- q 2 ) , 

•^bi-gi|>a- 1/4 

where we defined 

01,a = 2 f \V(x)\, 
J\x\>a-V* 

and we bounded the first integral analogously as we did with the integral in (jlO.lOp . Note that 
Pi, a — ► as a — > 0, because V G L 1 (M 3 ). We still need to control the last integral, on the r.h.s. of 
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the last equation. To this end, we observe that 

/ dpidp 2 dqidq 2 \<p(pi,P2)\ hKgi,92)| <S(pi + £>2 - 9i - 92) 

J|pi-9l|>Q- 1/4 

< 2 / dpidp 2 dqidq 2 \<p(pi,P2)\ |VKgi,92)| Hpi + P2 - 9i - 92) 

J\q 1 \>a- 1 / 4 /8 

f ^ (10-13) 

+ / dpidp 2 dqidq 2 \cp(pi,p 2 )\\ip(qi,q 2 )\5(pi + p 2 - qi - q 2 ) 

J\q 2 \> a -V*/% 



+ / | Pl |>a-V4 /2 dpid^dgid^l^Oi,^)! 1^(91,92)1^1 +P2 - - 92) • 

ki+92|<a- 1 /4/ 4 

The first two terms can be bounded by 

dp!dp 2 dqidq 2 \fi(pi,p 2 )\ 1^(91,92)1 &{Pi + P2 - 9i - 92) 

\qj\>a-V*/& 

< Ca 1 ' 12 {<p, ((Vi - V 2 ) 4 + (Vi + V 2 ) 2 + 1) ^> 1/2 ((Vi • V 2 ) 4 - Ax - A 2 + l) <^) 1/2 

(10.14) 

which holds for both j = 1, 2, and for a universal constant C, independent of a, cp, ip. To show (|10. 14|) 
note that, proceeding as in (|1U.1U|) (for example for j = 1), we have 



dpidp 2 d<?idg 2 \(p{pi,p 2 )\ 1^(91,92)1 5(pi + p 2 - qi - q 2 ) 
|gi|>a-V4/8 

< (<p, ((Vi - V 2 ) 4 + (Vi + V 2 ) 2 + 1) ^> 1/2 ((Vi • V 2 ) 4 - Ax - A 2 + 1) ^> 1/2 

/ \ 1/2 / \ 1/2 

dp \ I f dq \ 



X ' ; (p-(q-p)) 2 +p 2 + (q-p) 2 + l U&l>J\ q \> a -V*/8 \q~p\ A +p 2 + l 



(10.15) 

and thus (|10.14j) follows from (|10.4p and 

M , ,„ 1/4x1/3 



sup / 1 nj^— 5 < (8a 1/4 ) 1/3 sup /dg, 

P eR3 J\q\> a -i/y8 \q-p\ A +p 2 + l peR 37 |g-p| 4 +^ + l 



/ [g + pl 1 / 3 \ /• dg (10.16) 

V,SS» (|9| 4 +P 2 + 1) 1/6 J 7 (|9| 4 + 1) 5 / 6 



< (8a 1 / 4 ) 1 / 3 

< Ca 1 / 12 . 

As for the last term on the r.h.s. of (|10.13p . we note that 

/ | Pl |> Q -i/4 /2 d Pi d P2dqidq 2 \(p(p 1 ,p 2 )\\ ; tp(q 1 ,q 2 )\5(pi+p 2 -qi-q 2 ) 

< (<p, ((Vi - V 2 ) 4 + (Vi + V 2 ) 2 + 1) ^> 1/2 ((Vi • V 2 ) 4 - Ai - A 2 + 1) (^> 1/2 (10.17) 

/ r dq \ 1/2 f r dp \ 1/2 



\ q _ p |4 + p 2 + ! J I - y |p| »^i /2 (p . (5 - p))2 + p2 + (? _ p) 2 + 1 
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Since 



/ W>2| 



dp 

lpl - 2kl (p-(q-p)) 2 +p 2 + (q-p) 2 + l 



'|p|>a-i/4/2 ! 

| p |l/3 



< (2a 1 / 4 ) 1 / 3 I dp 
J\v\>2\a\ 



\p\>M (( p -iy-al\ +p 2 + ( q _ p) 2 + 1 



< Ca ^i dp, kp£H ( 10 - 18 ) 



<Ca^fd PjL ^ 



>M (p 2 -3l) +pl + 3± + l 
\P\ 1/2 



\p\ A + 1 
< Ca 1 ' 12 , 

it follows that the last term on the r.h.s. of (I10.13P is bounded by 

| P1 l> Q -i/4 /2 dpidp 2 dqidq 2 |£(pi,p 2 )| |?(?l,?2)| <*(Pl + P2~ Ql ~ 92) 
|gi+g 2 |<a-V4/2 

< Ca 1 / 12 (p, ((Vi - V 2 ) 4 + (Vi + V 2 ) 2 + 1) v9> 1/2 (V, ((Vi • V 2 ) 4 - Ai - A 2 + l) <^) 1/2 . 

(10.19) 

From the last equation, (|10.14[> . (jl0.13j> . and (I10.12j) . it follows that 

I II| < C(/3i )Q + a 1 / 12 ) 

x ((Vi - V 2 ) 4 + (Vi + V 2 ) 2 + 1) 99> 1/2 ((Vi • V 2 ) 4 — Ax — A 2 + l) V9> 1/2 • 



This together with f|10.11j) . implies ([HTH]I with /3 Q = C(/?i, a + a 1 / 12 + a 1 / 2 ). □ 

When dealing with the limit points 7^ t , for which we have stronger a-priori estimates, we will 
make use of the following lemma, whose proof can be found in [11] (Lemma 8.2). 



Lemma 10.3. Suppose that S a (x) is a function satisfying < S a (x) < Ca _3 l(|a;| < a) and 
J S a (x)dx = 1 (for example 8 a {x) = a~ 3 g(x/a), for a bounded probability density g(x) supported 
in {x : \x\ < 1} ). Moreover, for £ fC k , and for j = 1, . . . , k, we define the norm 

I Ji% := sup (X!) 4 . . . <z fc ) Vi> 4 . . . <*' fc > 4 (| jW(x fc ;x' fc )| + |V„ J (fc) (x fe ;4)| + |V^ j( fc )(x fc ; x' fc )|) 



k 



(10.20) 

and Sj = (1 — A x .) (here (x) 2 := 1 + x 2 ). Then if j( k+1 ^ (~x. k+ i,x.' k+1 ) is the kernel of a density matrix 
on L 2 (R 3 ( fc+1 )) ; we have, for any j < k, 

J dx fc+ idx' fc+1 J (fc) (x fc ;x' fc ) (5 ai (x k+ i - x' k+1 )5 a2 (xj - x k +i) - 6(x k +i - x k+1 )5(xj - x k+1 )) 

x 7 ( fe+1 )(x fc+1 ;x' fc+1 ) 

< C k \\J^l ( ai + v^) Tr|5 J 5 fc+l7 (fc+1) ^5 fe+1 | . (10.21) 
The same bound holds if Xj is replaced with x'j in 1110.21]) by symmetry. 
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A Properties of the cutoff function 0. 

Recall that the cutoff functions ejf ) = 9^ n) (x) defined for k = 1, . . . ,N and n € N, in Eq. fT4l . 
In the following lemma, whose proof can be found in [12\ Appendix A], we collect some of their 
important properties which were used in the energy estimate, Proposition 17.11 

Lemma A.l. i) The functions @^ are monotonic in both indices, that is for any n, k £ N, 

ejw < ei n) < i , e[ n+1) < < i . 

(n) 

Moreover, Q k is permutation symmetric in the first k and the last N — k variables, 
ii) For every k = 1, . . . , N , n S N we have 

n |v-e (n) 



3=1 u k 



Hi) For every fixed k = 1, . . . , N and n £ N we have 



(n) 
k 



< cr z e 



(A.2) 



B Removal of the assumption on derivatives of V 

The goal of this appendix is to explain how the assumption 

\V a V(x)\ < C for all x £ M 3 , \a\ < 2 (B.3) 

in Theorem 13,11 can be removed. The main observation is that (|B.3h is only used in the proof of 
the higher order energy estimate, Proposition 17. l\ in the form ||VViv||oo < CN 3 , ||V 2 VAf||oo < CiV 4 . 
More precisely, the estimate on ||VVaH|oo i s ^ TS ^ used in the study of the third term on the r.h.s. 
of (|7.19p (the third term on the r.h.s. of (|7.7p in the case k = 2); namely the term containing the 
commutator [T 1//2 ,L>fe] = [(-Hjv + A r ) 1 / 2 , Vi . . ■ Vfc]. Bounds on the first and second derivatives are 
also used in the proof of Lemma 17.21 However, in both cases, the final estimates turn out to be 
sub exponentially small in N (see (I7.23P and (|7.27p ). For this reason, the proof of Proposition 17.11 
remains unchanged if, instead of (|B.3p . we allow V = V^ N ' to depend on ./V and only assume 

||V a yW|U <e cN \ |a|<2, (B.4) 

for some sufficiently small k > 0. 

More precisely, suppose that the potential V > satisfies V(x) < C(x)~ a for some a > 5, with 
no assumptions on the derivatives V a V, for \a\ > 1. Then consider the evolution ipN,t = e~ lHNt i/)N 
of an initial iV-body wave function i/jn satisfying the two assumptions (13. ip and (13. 2p . with respect 
to the evolution generated by the Hamiltonian 

N N 

H N = ^ -Aj + ^ V N (xi - Xj) 

3=1 i<3 
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with V N {x) = N 2 V(Nx). As in Theorem O we claim that, for every fixed t G R, and every k > 1, 

(k) 

the /c-particle marginal , y N t associated with ipN,t is such that 

i%l^\n){vt\ m (B.5) 

as N — ► oo, w.r.t. to the trace-norm topology. 

To prove (|B.5p we can assume, without loss of generality, that the initial data ipjy is such that 

& n ,HUn) <C k N k . (B.6) 

In fact, if this is not the case, we can use the argument outlined in Section 2] (in the proof of Theorem 
13. 1 f) and based on the analysis of Section [9] (which does not use any assumption on the derivatives 
of V) to approximate tp^. 

From Theorem 16.11 it follows that the sequence T^\ = {7^ t }fe>i is compact with respect to the 
product topology r pro d defined in Section [H If we could prove, similarly to Theorem 17.31 that an 

arbitrary limit point = {7o^\}fc>i satisfies the a-priori estimates (I7.28p . it would follow from 
Theorem O that 

Too,* is a solution to the infinite hierarchy (|4.4p and, by the uniqueness result of 
Theorem 14.11 we could conclude the proof of (|B.5h using the same strategy outlined in Section [H 

To prove that every limit point t satisfies the a-priori estimates (|7.28p . we introduce a potential 
y( N ) = V*v( N \ where u^ N \x) = (const) e 3NK / 2 exp(— e N * x 2 ) with some sufficiently small k > (here 
the constant is chosen so that J dxv^ N \x) = 1), and we consider the evolution ipN,t = e~ lHNt \(jN of 
the initial data ^tv with respect to the modified Hamiltonian 

N N 
j=l i<j 

The potential 7^ satisfies the bounds \\V a V^ [| £ » < Ce^ m ^ 2 for all \a\ < 2. As we remarked 
above, this very weak control on the L°° norm is enough to prove Proposition 17.11 

Therefore, it follows from Theorem 17.31 that for any fixed t £ 1 every limit point T^t = {7oo\}fc>i 
of the sequence Tnj = {7wt}jfeLi (w.r.t. to the product of the weak* topologies) satisfies the bound 

Tr (1 - Ax) ... (1 - A k )j£] t < C k (B.7) 

for all k > 1. To show that a limit point T^j = {7^\}fc>i of the original sequence T^,t = {7^|}fc>i 
also satisfies this bound, it is therefore enough to prove that, for every fixed t E R, 

l^-felh ( B - 8 ) 

as N — > oo; in fact (1B.8I) immediately implies that every limit point 7^\ of Jm* is also a limit point 
of the sequence and therefore satisfies ()B.7p . 
To verify (|B.8j) . we observe that 
d 



dt 



ipN,t ~ i>N,t = 2 Im (ip N ,t, (Hn ~ H N ^j ipN,t) 

= N 3 (N - 1) Im (Vw,i, {V*(6 — v N )) {N(xi - x 2 ))$N,t) 

= N 3 (N - 1) Im I dxdy V^ t (x) V{y){5 - u^)(N( Xl - x 2 ) - i#jv,t(x) 



(B.9) 



(N - 1) Im / dy V(y) / dx iJ N Ax) (5 - v N ){x x - x 2 - y/N)^ N ^) 
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where we defined vn( x ) = N s u^ N ^ (Nx) (this implies that vn{x) = (const) N 3 e 3NK / 2 e n2 eN x2 is 
still normalized with H^ivHi = 1)- Therefore 



(B.10) 



(B.ll) 



< CiV||y||i sup (i>N,t, {$y - VN,y) (xi - x 2 )ipN,t) 
j/eM 3 

with 6 y (x) = 8{x — y) and v^,y{x) = vn{x — y)- It is simple to check that 

SUp (lpN,t, (Sy ~ VN,y) {x\ ~ X 2 )^N,t) 
y£R 3 

< Ce-TT(^ it , (1 - Ai)(l - A 2 )il> Njt ) 1/2 $N,t, (1 - Ai)(l - A 2 )^v, t > 1/2 . 
In fact, (jB.llh can be proven using the Fourier representation 

= J dp N - 2 dpidp 2 d qi dq 2 Sip! + P2 -qi- q^e^^-n) (l - e'^^ (Pi"9i) 2 /4 

V(^ + i)(p| + i) - , ^( q { + i)( q j + i) j 

with pat_2 = (P3, ■ ■ ■ ,Pn)- Using that |1 — e~ a \ < Ca 1 '^ for a > 0, applying a Schwarz inequality, 
and changing variables q — » p in V'TVt; we obtain 



(B.12) 



sup 



{lpN,t, (S y - VN,y){x\ ~ X 2 )i>N,t) 



< Ce s / dpAr_ 2 dpidp 2 dgidg2 S(pi + P2 - Qi - 92)- 



ki| 1/2 + l 



(qf + l)fe 2 + 1) 

X (Pi + l)(p| + 1) ( /?l^7V,t(pi,P2,PiV-2)| 2 + (3~ 1 $ Nit (pi,p 2 ,PN-2)\ 2 



for every f3 > 0. This implies (jB.lip because 



/ , |g| 1/2 + l 



< 00 . 



(B.13) 



(B.14) 



To bound the expectations of (1 — Ai)(l — A 2 ) on the r.h.s. of (|B,lip we observe that, as an 
operator inequality on L^(lR 3Ar ), we have 



N 



(H N + N) 2 = £(1 - Aj) + VN(xi - Xj) 



J =1 



2 



y=i / \i<j 

> N(N - 1)(1 - Ai)(l - A 2 ) - iV 6 



(B.15) 



X{ Xj ) 
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which, by (|B.6p . implies that 

Gtor.t, (1 - Ai)(l - A 2 )^, t ) < (iP N>t , ( ^NiN^l) ^ ) te ) 



(B.16) 



/ (gy + AQ 2 + CW 8 \ 
V N(N-l) J 



Using (|B.6P and a Schwarz inequality similar to (|B.15|) to compare and H™, it is simple to check 
that 

(i> N ,H%i> N ) <CN S 
and therefore, proceeding analogously to (|B.15P and (|B.16p . we also obtain that 

$ N ,t, (1 - Ai)(l - A 2 )^r,t> < CN 6 . (B.17) 
Inserting (|B.16[) and (|B.17[) into (jB.lip . and using (jB.lOp . we find 



d_ 

di 



ipN.t ~ 4>N,t\\ 2 



for some c > 0, which implies that 



1pN,t ~ IpN,, 



< Ct l / 2 e~ cNK (B.18) 



as iV — > oo, for every fixed t £ R. This completes the proof of (|B.8p . 

References 

[1] Adami, R.; Bardos, C; Golse, F.; Teta, A.: Towards a rigorous derivation of the cubic nonlinear 
Schrodinger equation in dimension one. Asymptot. Anal. 40 (2004), no. 2, 93-108. 

[2] Adami, R.; Golse, F.; Teta, A.: Rigorous derivation of the cubic NLS in dimension one. J. Stat. 
Phys. 127 (2007), no. 6, 1193-1220. 

[3] Anderson, M.H.; Ensher, J.R.; Matthews, M.R.; Wieman, C.E.; Cornell, E.A.; Science (269), 
198 (1995). 

[4] Bardos, C; Golse, F.; Mauser, N.: Weak coupling limit of the A-particle Schrodinger equation. 
Methods Appl. Anal. 7 (2000), 275-293. 

[5] Bourgain, J.: Global well-posednedd of defocusing 3D critical NLS in the radial case, Jour. 
Amer. Math. Soc. 12 (1999), 145-171. 

[6] Colliander, J.; Keel, M,; Staffilani, G.; Takaoka, H., Tao, T.: Global well-posedness and scatter- 
ing for the energy-critical nonlinear Schrodinger equation in M 3 . Ann. of Math. (2) 167 (2008), 
no. 3, 767-865. 

[7] Davis, K.B.; Mewes, M.-O.; Andrews, M. R.; van Druten, N. J.; Durfee, D. S.; Kurn, D. M.; 
Ketterle, W.; Phys. Rev. Lett. (75), 3969 (1995). 

[8] Elgart, A.; Erdos, L.; Schlein, B.; Yau, H.-T. Gross-Pitaevskii equation as the mean filed limit 
of weakly coupled bosons. Arch. Rat. Mech. Anal. 179 (2006), no. 2, 265-283. 



51 



[9] Elgart, A.; Schlein, B.: Mean Field Dynamics of Boson Stars. Commun. Pure Appl. Math. 60 
(2007), no. 4, 500-545. 

[10] Erdos, L.; Schlein, B.; Yau, H.-T.: Derivation of the Gross-Pitaevskii Hierarchy for the Dynam- 
ics of Bose-Einstein Condensate. Commun. Pure Appl. Math. 59 (2006), no. 12, 1659-1741. 

[11] Erdos, L.; Schlein, B.; Yau, H.-T.: Derivation of the cubic non-linear Schrodinger equation from 
quantum dynamics of many-body systems. Invent. Math. 167 (2007), 515-614. 

[12] Erdos, L.; Schlein, B.; Yau, H.-T.: Derivation of the Gross-Pitaevskii Equation for the Dynamics 



of Bose-Einstein Condensate. Preprint arXiv:math-ph/0606017. To appear in Ann. of Math. 



[13] Erdos, L.; Yau, H.-T.: Derivation of the nonlinear Schrodinger equation from a many body 
Coulomb system. Adv. Theor. Math. Phys. 5 (2001), no. 6, 1169-1205. 

[14] Ginibre, J.; Velo, G.: The classical field limit of scattering theory for nonrelativistic many-boson 
systems. I.-II. Comm. Math. Phys. 66 (1979), no. 1, 37-76, and 68 (1979), no. 1, 45-68. 

[15] Ginibre, J.; Velo, G.: On a class of nonlinear Schrodinger equations with nonlocal interactions. 
Math Z. 170 (1980), 109-136. 

[16] Ginibre, J.; Velo, G.: Scattering theory in the energy space for a class of nonlinear Schrodinger 
equations. J. Math. Pures Appl. 64 (1985), 363-401. 

[17] Hepp, K.: The classical limit for quantum mechanical correlation functions. Comm. Math. Phys. 
35 (1974), 265-277. 

[18] Lieb, E.H.; Seiringer, R.: Proof of Bose-Einstein condensation for dilute trapped gases. Phys. 
Rev. Lett. 88 (2002), 170409-1-4. 

[19] Lieb, E.H.; Seiringer, R.; Solovej, J. P.; Yngvason, J.: The mathematics of the Bose gas and its 
condensation. Oberwolfach Seminars, 34. Birkhauser Verlag, Basel, 2005. 

[20] Lieb, E.H.; Seiringer, R.; Yngvason, J.: Bosons in a trap: a rigorous derivation of the Gross- 
Pitaevskii energy functional. Phys. Rev A 61 (2000), 043602. 

[21] Kato, T.: On nonlinear Schrodinger equations, Ann. Inst. H. Poincare Phys. Theor. 46 (1987), 
113-129. 

[22] Klainerman, S.; Machedon, M.: On the uniqueness of solutions to the Gross-Pitaevskii hierarchy. 
Commun. Math. Phys. 279 (2008), no. 1, 169-185. 

[23] Michelangeli, A.: Equivalent definitions of asymptotic 100% BEC. // Nuovo Cimento B 123 
(2008), no. 2, 181-192. 

[24] Reed, M.; Simon, B.: Methods of modern mathematical physics: Scattering Theory. Volume 3. 
Academic Press, 1979. 

[25] Rodnianski, I.; Schlein, B.: Quantum fluctuations and rate of convergence towards mean field 
dynamics. Preprint arXiv:math-ph/071 1.3087. To appear in Comm. Math. Phys. 

[26] Rudin, W.: Functional analysis. McGraw-Hill Series in Higher Mathematics, McGraw-Hill 
Book Co., New York, 1973. 



52 



[27] Spohn, H.: Kinetic Equations from Hamiltonian Dynamics. Rev. Mod. Phys. 52 (1980), no. 3, 
569-615. 

[28] Strauss, W.: Non-linear scattering theory at low energy. J. Fund. Anal. 41 (1981), 110-133. 

[29] Yajima, K.: The H^ fc ' p -continuity of wave operators for Schrodinger operators. J. Math. Soc. 
Japan 47 (1995), no. 3, 551-581. 

[30] Yajima, K.: The H^ fc ' p -continuity of wave operators for Schrodinger operators. Proc. Japan 
Acad. Ser. A Math. Sci. 69 (1993), no. 4, 94-98. 



53 



